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DEVELOPMENT OF A NONLINEAR VORTEX METHOD 
By 


Osama A. Kandil* 


ILis semi-annual report covers the progress of the research *«rk 
conducted under this grant from October L, I9H0 to March 31, 1981. IXiring 
this period, the following tasks have been accomplished: 

1. Development of the Modified Nonlinear Discrete Vortex (MNDV) method 
has been completed. Ihe method is complemented by a viscous model for the 
vortex core based on the first-order boundary-layer-like equations. (Viscous 
analysis is supported under a separate contract from the Naval Air Develop- 
ment Center, Warminster, Pa.). Tlie results of the inviscid and viscous 
methods have been reported in the AIAA paper No. 81-1263 titled "Recent 
Improvements in the Prediction of the Leading and Trailing Edge Vortex Cores 
of Delta Wings." A copy of this paper is enclosed in Appendix A. This 
paper has been submitted for publication in the Journal of Aircraft. 

2. Development of the Nonlinear Hybrid-Vortex (NHV) method for wings 
with side-edge separation in steady and unsteady flows has been completed. 
The steady and unsteady results of this method have been reported in the 
AIAA paper No. 82-0351 titled "Steady and Unsteady Nonlinear Hybrid Vortex 
Method for Lifting Surfaces at Large Angles of Attack." A copy of this 
paper is also enclosed in Appendix A. This paper has been submitted for 
publication in the AIAA Journal. 

3. The principal investigator has been invited as a lecturer in a 
short course on "Computational Methods in Potential Aerodynamics," organized 
by the International Center for Transportation Studies (ICTS) in Amalfi, 
Italy, in the period from May 31 to June 5, 1982. The lecture is titled 
"Steady and Unsteady Incompressible Free-Wake Analysis." The lecture covers 
the state of the art of free-wake analysis including a detailed presentation 


^Associate Professor, Department of Mechanical Engineering and Mechanics, 
Old Dominion University, Norfolk, Virginia 23508. 


of Che work developed so far under this granc. A copy of Che lecture is 
also included in Appendix A. The lecture will appear in the proceedings of 
Che course, an ICTS publication. 

The following Casks are to be completed by September 30, 1982: 

1. The Inverse Hybrid-Vortex Method; Accurate Near-Field Calculations. 
For accurate near-field calculations in the MNDV method, the lumped vortex 
filaments, representing the wing and its wake, are redistributed on Che 
corresponding surface panels. This modification guarantees Che accuracy of 
Che MNDV method and moreover maintains the efficiency of this method over 
Che expensive panel methods. To eliminate any leapfrogging at the location 
of crailing-edge core and to expedite the convergence, functions of helical 
curves are used to fit the computed nodes of Che free-vorCex segments after 
each iteration cycle. 

2. Steady and Unsteady Nonlinear Hybrid-Vortex Method for Leading-Edge 
Separation. The present NHV-method (covered in AIAA Paper Mo. 82-0351) is 
currently extended to treat wings with leading-edge separation. Ihe 
computer code is very general; it can treat combined-edge separation (LE & 
TE), general planforms, and cambered wings. The solenoidal property of 
vorticity is also enforced across the edges of vortex panels. This 
eliminates (on the average, since the vorticity distribution is assumed to 
be linear) the discontinuity in the vorticity gradient. 

3. The steady part of the NHV-method completed in item 2 will be used 
to develop the unsteady NHV-method for the frequency domain approach. 

4. The NHV-method for bodies of revolution with prescribed separation 
lines is being completed. 

Mr. Li-Chuan Chu is working on the NHV-method reported in items 2 and 3 
above for his Ph.D. dissertation while Mr. Chong is working on the NHV- 
method reported in item 4. 
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/APPENDIX A 

1. AIAA Paper No. 81-1263. 

AIAA Paper No. 82-0351. 

3. Lecture preuented at the International School of Applied 
Aerodynamics, Amalfi, Italy, Hay 29-Jun-* 3, 1982. 
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STrjkOY AND UNSTEADY NONLINEAR HYBRID VORTEX NETHOU 
FOR LIFTING SURFACES AT LARGE ANGLES OF AHACK^ 

OtMM A. Kandll*, Ll-Chu«n Chu** and Thonas Turuud 
Old Dominion Unlvtrslty, Norfolk, VA. 23S0B 


Abstract 

Sttady and unstaady !(on1 Intar Hybrid Vorttx 
(NHV} ntthod, for 1 om asptct ratio wings at largt 
angits of attack, Is dtvoloptd. Tht ntthod osts 
vorttx ptntls with f1rst*ordtr vorticity distribu* 
tion (tquivaltnt to stcond*ordtr doubitt distribu- 
tion) to calculatt tht Inductd vtlocity In tht nttr 
fltld using clostd fona txprtsslons. In tht far 
fitid, tht distrlbuttd vorticity Is rtductd to con- 
ctntrattd vorttx lints and tht slmpltr Blot-Savart's 
law Is anploytd. Tht sMthod Is appiltd to rtctan- 
gular wings In sttady and unstaady flows without 
any rtstrictlon on tht ordar of magnitudt of tht 
dlsturbancas In tht flow fltld. Tht numtrical 
rtsults show that tht ntthod accurataly prtdicts 
tht distrlbuttd atrodynaailc loads and that It Is of 
acctptablt computational tfficitncy. 


I. Introduction and Background 

In rtctnt ytars, dtvtiopmtnt of numtrical 
mtthods for prtdicting tht sttady and unstaady 
atrodynamic characttrl sties of lifting surfacts 
txhibiting loading- and/or sidt-tdgt stparatlons 
has rtctivtd considtrabit atttntlon. 

' For tht sttady-flow probitms stvtral numtrical 
ttchniquts havt bttn devtloptd. Thtst Includt tht 
NonI Intar 01 serttt- Vorttx^ (NOV) Mthods, hlgh- 
ordor doubitt panti mtthoUs'~* and Nonllnaar Hybrld- 
Vorttx CMHV) mtthods. '0-11 For tht unsttady-flow 
probitms, tht lltoraturt shows ftwtr rEumtrlcal 
ttchniquts which Includt tht NDY-mathodslZ-IS and 
constant doubitt panti nothodslO-17, Tht littra- 
turt lacks high-ordtr panti mtthods for tht unstaady 
^ow problems. 

For this reason and because of the success of 
the high-ordtr panel mtthods In sttady-flow pro- 
blems, wt art presenting In this paper an efficient 
and accurate method for tht steady- and unsteady- 
flow probitms of lifting surfacts at large angits 
of attack. 

In this method, vorttx panels with first-ordtr 
vorticity distribution Is used In the ntar-fleld 
calculations. In tht far field calculations, the 
distributed vorticity over each far- field panel Is 
lumped Into equivalent concentrated vorttx lines. 

In this way, accuracy Is satisfied In the near 
field while computational efficiency Is maintained 
In tht far field. Tht coupling of a continuous 
vortex-sheet representation and a concentrated 
vorttx-llne representation for solving tht non- 
linear Hftlng-surfact problem Is called the 


t This research work Is supported by NASA Langley 
Research Center under Grant No. NSG 1560 and by 
the Naval Air Development Center under Contract 
No. N62269-80-C-0704. 

* Associate Professor of Mechanical Engineering and 
hicchanics, AIAA Kambcr 
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"Nonlinear Hybrid Vortex (NHV)" mettod. 


. II. Formulation of the Problem 

The problem Is formulated relative to a wing- 
fixed frame of reference xyz. The x-axIs Is the 
wing centerline and the xy-plane Is the wing plene 
of syamMtry. The wing is rotating at the angular 
velocity 0 and the freestream velocity Is u . The 
general orientation of the wing Is described using 
the Eulerlan angles a, B, and y, figure 1, which 
refer to the angles of attack, yaw and rollj respec- 
tively. In terms of these angles and their time 
rete of change, the dimensionless freestream velo- 
city and the angular velocity are expressed by 

7 • cosa cpsB T ♦ (-slna cosy COSO SinB s1ny)J 
" ♦ (slna siny + cosa sing cosy)k (1) 


0 ■ (;a sins ♦ y)T ♦ (a cosg siny ♦ S cosy)J 
+ (a COSB cosy - B Slny)k ■ ♦ 0^,7 (2) 

where T, J and 7 are the base unit vectors of the 
xyz- frame of reference . 

The unsteady Irrotatlonal Ideal flow In thr 
region outside of the wing and Its separated vortex 
sheets Is governed by the Laplace's equation 

■ 0 (3) 

where a(7,t) Is the perturbation velocity potential. 
The no-penetration condition on the wing s(f,t), 
relative to the wing fixed frame of reference. Is 
given by 

♦ (7_ + 7a - ax7)*7$ *0 on s(7,t) ■ O (4-a) 

For a rigid wing, |f - 0 and equation (4-a) reduces 
to " 

(7^ + 7a - SxF)*7j *0 on s(7) • 0 (4-b) 

On the separated free-vortex sheet w(?,t), the no- 
penetration condition Is given by 

♦ (7^ ♦ 7a - ox7)*iw • 0 on w(7,t) ■ 0 (5) 

The no-pressure Jump condition on w(7,t) Is obtained 
from the unsteady Bemoulll's equation 

Cp(7,t) ■ -7a*[Y4 ♦ 2(7. - 0x7)] - 2 H (a) 

where ^(7,t) Is the pressure coefficient at any 
point rand at any tine t. Forming the pressure 
Jump from equation (a) and equating the result to 
zero, we obtain 

dtp • Cp^ - Cp2 - -(7a^ - 7a2)-[7a,+7a2*z(V«^^n 

. 2 |i^i * ♦zJ . 0 (b) 
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OF POOR QUALITY 


MtMTt tht subscripts 1 and Z rsftr to tho uppor and 
loiMr surfacts of tho wing; rospoctlvaly. Roarran- 
glng aquation (b) and totting 

a^ - #2 ■ ^ (c) 

ono obtains 

ACp - -2C|^ ♦ "7^ • 7)Ua) ■ -2 - 0 (d) 

whoro 7 ■ volocity of a wofco olaaont rolativo to 
tho xyx'fraao of rofortnco 

• j * 7*2^ ♦ • S*F 

Eguotlon (d) roprttonts tbo tboorou of Klolvin and 
Holariioltz 

^ ■ ^ // w • n^ d A • 0 on w(7,t) ■ 0 (6) 

of consorvotlon sf tbo circulation and tbo outflow 
of vorticity; rospoctivoly. In aguotlon (6), n. It 
a unit nor«a1 to tbo surfaco A boundod by a cloiod 
curvo around whieb tho circulation r It calculatod. 
Equation (6} slaply statos that tbo rate of chango 
of circulation around a closad curvo or tho rato 
of chango of outflow of vorticity tbriugb tho sur- 
faco boundod by this closad curvo Is zoro (following 
tho sasM fluid particlosj 

For uniquonoss of tho solution, ono bat to 
Iwposo tho Kutta condition along tho adgos of 
soparatlon. Haro, Kutta condition Is roprosantad 
by 



Finally, tha Infinity condition roquiros that 

0 away froa s and w (8) 

Equation (1} - (8) arc tho roquirod aquations for 
tho gonoral unstoady flow problow. 


III. NonI Inoar Hvbrid-Vortos Hothod for Stoadv Flows 
III.I. Govomlnq Equations 


Tho rosoarch offort In this papar Is concon- 
tratad on tha sjosMtrlc-flow problaws. FOr a 
stoady-syanatrlc flow, tho govoming aquations aro 
obta1nad_frow tho gonoral aquation (I) • (8) by 
sotting Q-0, a-Y-Oand dropping tho tlao 
dopandont torn. Tho rosulting aquations an 


a, • cos a T - sin a 7 

• 0 

Ca, ♦ *a) • • 0 on s(r) • 0 

Cn, ♦ da) • 7^ • 0 on w(7) • 0 

ACp ■ •2(n’^jto) • (o^ ♦ da)" 0 on w(?) • 0 


‘SI ■ " 

Tt.LE.SE 

va * 0 away froa s :nd w 


( 9 ) 

( 10 ) 

(11) 

( 12 ) 

(13) 

( 14 ) 

(15) 


III.2. Mtthod of Solution of tha Stoadv Flow Froblaw 

Equation (12) roquiros tho flow to bo tangont 
to w whilo option (13) roquiros this tantontlal 
flow to bo parallol to tho vorticity direction. 
Thoroforo, If tho flow direction Is forced to bo 
parallol to tho vorticity dIrKtIon on the surfaco 
w, tho boundary conditions of aquations (12) and 

(13) aro autoasitlcally Mtisflod. Next, wo outline 
tha aathod of solution. 

Once tha wing and Its froa*shaar layers are 
roprasantad by vortex sbaats, aquations (10) and 
(18) are autaoMtlcally satisfied. Tho basic un- 
knowns *1n tha„prasant problon are tha vorticity 
distribution « and tha fraa-fortax shaat w. Thoy 
are datamlnad by satisfying tha rosHlnlng boundary 
conditions, aquations <11M 41, through a flnlta- 
olonont ty^ approach. 

In this approach, tha bound-vortax shaat 
(raprasanting tha wing) Is divided Into quadrilataral 
bound-vortex panels while tha fraa-vortax shaats 
(raprasantlno the froa-shaar layers) are divided 
Into triangular fraa-vortax panels, soa figure 2. 

On each vortex panel (bound or froa), a local 
vorticity distribution with undotomlnad coeffi- 
cients Is defined In a I oca I -coordinate systan. 

Tha local distribution Is saloctad such that tha 
solonoldal property of vorticity Is satisfied. Tho 
contliMity of vorticity (a coapatiblllty condition) 

Is anforcad at certain nodes on tho Intaroloaort 
boundaries of adjacent panels. 

Tha roaolning boundary conditions, aquations 
OV-04), are enforced at certain nodes of t«o vortex 
panels to obtain tha undetarwlnad coafflcants of 
tha local vorticity distribution and tha shape of 
tha fr aa -vortax panels. Kutta condition, aquation 

(14) , Is onfOrcad at the nodes of the bound-vertex 
panels along tho edges of separation. Tha no- 
panotratlon condition, option (II), Is onfOrcad 
at tho average points of tho bound-vortex panels. 

Tho no-ponotratlon and no-prossuro-Juaip conditions, 
equations (12) and (13), aro slaultanoously satis- 
fied at tbo nodes of tho froa-vortox panels. 

To satisfy those conditions, an Itarativo 
tochniqua Is followed which altomativoly yields 
tha local vorticity distribution on tho bound- 
vortex panels and tho shape of tho freo-vortox 
panels. During a typical Iterative cycle, an 
ovordotoralnod set of algebraic equations aro solved 
for tho undotonalnod coefficients of tho local 
vorticity distributions. This Is followed by adjus- 
ting tha surface w such that v and (T_ ^ Vb) at tho 
nodes of tbo freo-vortox panels bacoaa parallol. 

Tha ovordatomlnad tot of aquations consists of 
tha equations obtained frtxa tha contlmilty of 
vorticity condition, Kutta condition, tha no- 
panatratlon condition on tha wing, and a syMMtry 
condition appHod at tha nodes along tha Una of 
tymatry. This sot It solved by a Taast-squaro 
tochniqua. 

Once tho Iterative technique converges, tho 
pressure distribution Is calculated and this Is 
followed by calculating tho total aerodynaalc 
characteristics. Convergence of the technique Is 
expedited by using an Initial guess for w provided 
by tho Nonlinear 01 scroto-Vortox aethod4~4 (nOV- 
nothod). 
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In tht ntxt substctlons. tht basic aquations 
ptquiraq at tach stop of tht solution ara givan. 


1 1 1.2.1. Quadrilataral Panel 

Quadrilataral vortax pantls ara usad to repra- 
sant the bouno-yortex sheet. On each panel, a local 
first-order vorticity distribution is specified. 

For the KUl P^nel , the vorticity distribution Is 
given by 

■ •i(« ♦ dgOo c(io ♦ djOy c(w ( 16 ) 

-^(lO -*2(« C(« - e^OO - ajdO cUO (17) 

which contain five undetenelned coefficients; a.-ie* 
It should be noted that a first order vorticity' 
distribution Is equivalent to a second-or^er doublet 
distribution. The distributions given by equations 
(16) and (17) satisfy the solenoldal property of « 

7 • w ■ 0 (IB) 


OWGlNft'- 

OF POOR QUALITY 


where 


• (w^ sin tt 2 - CO* Y 2 )* 2 /(cos * 1 " ^2 

- sin cos Y 2 ) (21 ) 

■ (ij cos Y] - *1n Y^)t^/(cos Y^ *1n T 2 

- sin Y^ cos Y 2 ) 


,3 'i«-1 c(c) 

i ■ T I 7 ; w, dC <U 

' ^ 1 -1 c, 0 ' 


, 3 'l^l C(C) 


^ '^i-i c 


r ; f w, dc <u 


( 22 ) 

(23) 

(24) 


and A Is the panel area. Figure 3.b. shows the 
different parasMters used In equations ( 20 ) - ( 22 ). 


The four corners of the panel serve as nodal points 
where continuity of vorticity condition, Kutta con- 
dition, and/or syeseetry condition are satisfied. 

The c and ' axes are located In the psnel plane 
such that the ;- axis coincides with t'le 1-4 side 
of the panel. The n-axis Is perpendicular to the 
panel such that ;, t, and n axes fore a right-handed 
local coordinate system. The average point of the 
four nodes serves as the control point where the ^ 
no-penetration condition Is enforced, see figure 3a. 


In the neer-fleld, the Induced velocity at any 
field point X, y, z Is calculated by 

. ,2,3 H+1 c(C) y », e ♦[(z-c)w 

v(x.y.4) " it I 7 / .} - ^ ; 3/? 

^"1-1 0 [(x-c)SMz-c)^]^'^ 


-(*-e)wj7^ - yu. e, 
X — a 5 — *- (n dc 


(19) 


where the suoeatlon over 1 Indicates that the panel 
Is divided Into two (In the case of a triangular 
panel, see sec. II 1. 2. 2.^ or three subpanels. Sub- 
stituting equations (16) and (17) Into equation 
(19), we integrate the resulting expression In the 
(-direction analytically. The resulting equations 
are given in Appendix A. The Integration in the 
(-direction is perfomed numerically. Currently, 
complete analytical Integrations in both directions 
are developed. Appendix B. 


III. 2. 2. Triangular Panels 

Triangular panels model nonplanar and twisted 
surfaces more accurately than quadrilateral planar 
panels. Therefore, they are used to represent the 
free-vortex sheets where highly nonplanar and 
twisted surfaces are encountered. The local 
vorticity distributions are still giver, by equations 
(16) and (17). The corners of this panel serve as 
nodal points where certain boundary conditions Ire 
enforced. The c-ex1s coincides with the 1-3 side 
of the panel, the (- axis Is In the panel plane, 
and the n-ax 1 s Is perpendicular to the panel such 
that c, C, and n axes form a right-handed local 
coordinate system, figure 4a. 

In the near field, the Induced velocity at any 
point Is given by equation (19), where k ■ 1,2. In 
the far field, the vorticity distrloutlon on the 
panel Is lumped Into concentrated vortex lines so 
that 31ot-Savart' s law can be used to calculate the 
Induced velocity. The locations of the two concen- 
trated vortex lines are chosen to be the line con- 
necting the midpoints of sides 1-2 and 2-3 and the 
line connecting node 2 to the midpoint of side 1-3. 
The far-fleld equations of the Induced velocity are 
given by equations (20)-(24), where y, * 0 and 
k - 1 , 2 . ' 


I I I. 2. 3. Boundary Conditions 


In the far-fleld, the vorticity distribution 
Is 1 taped into concentrated vortex lines and the 
simpler Blot-Savart's law Is used to calculate the 
Induced velocity. The locations of the concen- 
trated vortex lines ere the lines connecting the 
mid-points of opposite sides of the panel. The 
far-fleld equation of the Induced velocity Is given 
hy 

V . ^ (cos e,. ♦ cos 62 j)e^^ ♦ (cos 

% ’ c 

♦ cos « 2 ()«y( (2C) 


At each control point of the bound-vortex 
panels, we satisfy equation (11). For the mtli 
control point (a receiver panel), equation (11) 
gives 

H _ 

r v(m,k) • 'n(m) - -e • n (m) (25) 

k-1 » - » 

where N is the total number Yf bound- and free- 
vortex panels (k refers to a sender panel], and 

5 

v(m,k) ■ r Cu,(m,k)e (k) ♦ v.(m,k)e (k) 

j.l J C J n 

♦ Wj(m,k)e,(k)]aj(m,k) (26) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Tlw coflfflclants of 0., 9., and w, art ^Ivan by tha 
aduotlont of Ap|iaiidlK''|A)''or AppaiMla (8). 

Tha unit noraal 17 (ii}jof a b 2 ynd-vortax panal 
and tlw basMunlt vdctort a.(k). • (k) and a (k) of 
tha local coordlnatai of tha k |h Bandar panal 
(bound or fraa) ara axprasiad In toms of tha bata* 
unit vactors of tha global Mlng-flxad coordinata 
lyttaa T. J fOlloM: 

nj(«) • 8 ^(«) T ♦ d^y(«) J ♦ ^ ( 27 ) 

•^(k) • dj^Ck) T ♦ dj^(k) T ♦ dj,(k) r 

Mliara d rafar to tha diractlen eealna . 


Nhara » and « ara axprattad In tanas of tha 
undartahalnad coafficlantt a(J,k) using aguatlons 
(16) and (17). 


Kutta condition, aq«iat1on (14), Is satisfied 
at tha nodes of bound*vortax panels along tha edges 
of separation. At a typical global node 02, J), 
aquation (14) gives 


•j,(l+2,J) v,(l*2,J) - -j(I+2,J) v^(Ia2.J) - 0 (31) 

Since m , V and v ara functions of tha indar- 
nlnad coaffTclaflts, equation (31) Is a ncaUrinar 
aquation In i(J,k). Adding aquations of this typo 
to aquations (29) and(31) disturbs tha linearity of 
tha resulting tat of aquations. Thorafora, wa 
anforca aquation (31) In a llnaarltad aporMch. 

Ua divide aquation (31) by v(I^2,J) and let 


*,(I+2,J)/v(I+2,J) - eosa(I+2,J) 


Substituting aquation (28) Into (26) and sub- 
stituting tha resulting aquation and aquation (27) 
Into aquation (2S), tha no-panatratlon at any nth 
point reduces to 


N S 
Z £ 
k«1 J«i 


-nx 

T 

K. 

“nx 

^c* 


“j 

“ny 


“cy 

'rvy 

“«y 



“nz 

t 

.“C2 



i 

> 


T 


cose 




nx 

-Sine 


d. 


ny 

0 




a 


( 29 ) 


It should bo noted hare that tha nunbar of 
unknowicoafficlants a(J,k) are reduced fron 5N 
(5x total nunbar of bound- and fraa-vortax panels) 
to S Nv (5x total nuober of bound-vortax panels) 
by equating tha outflow of vorticity fron tha bound- 
vortex panels to the Inflow of vorticity Into tha 
fraa-vortax panels along tha separation edges. 
Noraovar, tha vorticity vector Is aligned with one 
tide of tha triangular fraa-vortax panels. Once 
this Is acconpllshad, tha fraa trgunant n has tha 
range 1 to 8^. and thus wa obtain unknowns. 


In addition, we write aquations of vorticity 
continuity at each coaon node of tha bound-vortex 
panels. A typical couple of aquations written at 
the (Irl,J) olobal node, corresponding to the two 
local nodas (2,k) and (1,kr1) betw een panels nueber 
k and kri is given by (sea figure a) 

-jd^UJ) - -5(2,k)d^j(k) ♦ «j(2,k)d^j(k) 


• -;(l.krl)dj^(k+l) r »^(2.krl)dj^(k*1) 

(30 

-,(!♦! J) • -j(2.k)d^j,(k) ♦ «j(2,k)d^jj(k) 

• Wj(l ,krl )dj^{kr1 ) ♦ *^(1 ,krl )d5jj(krl ) 


and v^(I’^2,J)/v( 14>2,J) - Sln6(l-r2,j) (32) 

Thus, aquation (32) becoaes 

w,(I^2.J) C0S|(H2,J) - w,(H2,J) s1n6(I*2,J) -0 
* * (33) 

Assueing that 6(1*2, <)) Is Initially known, aquation 
(33) bacoaws a linear aquation Incas. In tha sub- 
sequent Iterative steps, the angle • Is calculated 
by aquations (32) using thaoBvaluas of tha pre- 
ceding Iterative step. 

Tha last set of linear equations ara obtained 
froe the syMotry condition (fpr syneatrlc flow) 
along tha. root chord. A typical svMotry condi- 
tion written at a global node (I,J), corraspendlng 
to a local node (1,£), Is given ^ 

-^(*I,J) • Wjd.kldjjjdl) ♦ •5(l,£)d^jj(B) - 0 (34) 

Tha resulting sat of aquations obtained free 
aquations (29). (30), (33) and (34) are salved for 
the coefficients ta) by a laast square solver (ovar- 
daterelnad sat of aquations). 

(hica tha coefficients ara datarelnad, aquations 
(12) and (13) ara enforced by aligning sides 1-3 of 
tha triangular fraa-vortax panels with tha local 
flow directions. Cations of the Allowing fore 
are used to calculate tha new global downstraae 
nodes of these sides: 

x(L*l,J) • x(L,J) ♦ &fgm(L,J)/Vj(L,J) 

y(L*l,J) - y(L,J) ♦ bv^(L,J)/v,(L,J) (36) 

t(L+l.J) • *(L,J) ♦ 5 v^(L,J)/v,(L,J) 

where (L,J) and (L*1,J) rafar to the upstraao and 
tewnstraao nodes of side (1-3), b Is tha length of 
this side, and v Is tha velocity calculated at tha 
upstraaia node (L,J). 


I I 1.2. 4. Pressure Distribution 

Tha nat-pressurw coefficient Is calculated at 
tha no-panatratlon control point of a bound-vortex 
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IMfttl. At tlM ■ th control point, tho not'prttsuro 
coefficient Is given by 

H N 

aC.(») ■ 2{w„(«) l V (e.k) - w,C»)C I V («,k) 

P * H.l * * l(.1 » 

* cos a]} (36) 

where «,,(«) • * Wj(«)djj,(e) 

wjt«) • WjOhldjjCe) ♦ 

Equations (If). (17), and (2f) give « . m and v 
In tenu of the coefficients a. The herall force 
coefficient Is than calculated by 

"b "b 

C • r aCpOi) A(u)/ t A(>) (37) 

where A(p) Is the area of • th panel. A no re 
accurate calculation of C Is obtained bv fitting 
a surface to discrete of Values of aC.(o). The 
resulting function aCp(x,z) Is then us^ to obtain 


C„ • ff aC_(x.z) dA//; dA (38) 

" A P A 

IY.1. Nonlinear Hybrid-Vortex Hethod for Unsteady 

TTfiS 


Pdr.uns(eady syuaetrlc flow, we set 
• ■Y*S*T"0 1n equations (1) and (2). The 
resulting equations are 


e^ • cos e T - sin o J 


(39) 

0 • Oj 7 ■ i r 


(40) 

v*e ■ 0 


(41) 

(e^ ♦ q> - Qjx7) • • 

0 on s(7)»0 

(42) 

If ♦ (a, ♦ fp - Oj*r) . 

Vw ■ 0 on w(7,t)“0 

(43) 


on w(7,t)a0 

(44) 

r> 


(45) 

TE,L£,SC 



tp - 0 

away froa s and w 

(46) 

IV. 2. Method of Solution of the Unsteady-Flow 


Problea 




The source of flow unsteadiness In this probloB 
can be a tiue dependent angle of attack or a tlue* 
dependent freestreoi speed. Figure t . In the tiue 
doaain approach of this orobleu, we divide the 
rHinctlon or angle of attack into discrete changes 
In the angle of attack corresponding to discrete 
changes In tlae; l.e., at t-tg, a-og and at t*tg 
* at, og^ etc. Tim preblan Is then solved at 
each tloe step where the solution of the precedir.g 
tlae step serves as the Initial condition for the 
present tlae step. 


This approach can treat problau where the flew 
unsteadiness starts (at t • t ) froai a steady flew 
or where the flow unsteadiness starts lapulslvely 
froai rest. The fonaer problaai Initially requires 
the steady-flM solution to be known; a^tiuns (9)> 
(15), while the latter problaai Initially f«<3>t<res 
tha solution of the flow over a wing without a 
wake surface. 

Once the Initial condition Is obtained, we 
•arch stop by step In tlae Mtisfying aquations 
14^) - (45). J^t each tlae step, the basic unknowns 
w(r,t) and w(r,t) are ebulned. 

Next, wa discuss the bewdary conditions and 
the calculations of the pressure distribution at 
each step. 


tLgili , 

Although the developed t1ne>*d03n1n technique Is 
a general one and ■oreover Is net restricted to any 
particular source of flow unsteadiness, the case of 
an lapulsively started wing frrai rest Is censIderetJ 
for the purpose of explaining tha,deta11s of the 
technlqjtjB. In this case, we set wra and replace 
T by -e_ In equations (39), (40) end (42)-U5). The 
following steps explain the procedure to enftree 
the boun^£y conditions of equations (42)>(4S) and 
to obtain w(7,t) and w(r,t): 

a. The Initial condition, at t«t , Is considered 
to be a wing without a wake surface having a velo- 
city of -e,. At this Instant, we assuae that net 
enough tlae has passed for the vortlclty to be 
convected froa the separation edges^ The Initial 
vortlclty distribution on the wing wCr.t) Is ob- 
tained froa the least-square solution of an ever- 
detenslnad tat of linear a I Mbra 1c actions In the 
unknown coeffIciMitt (a). The tat of equations 
consists of the no-penetraclon conditions, equation 
(42), the continuity of vortlclty Mnditlont, 
equat1on(30), the Kutta conditions, equation (45), 
and the syaaetry conditions, equation (34). 

It should be noted that the Initial vortlclty 
distribution on the wing Is such that the circula- 
tion around any closed curve Mbracing a wing sec- 
tion or equivalently tha outflow of vortlclty froa 
a surface enclosed by this curve Is xero. Con- 
sequently, a starting vortex of opposite strength 
to that of the vortlclty on the wing develops at 
the edges of separation. The starting vertex Is 
then convected downstreaa with the local particia 
velocity, the flew existing at the preceding Instant 
Is disturbed, and the vortlclty distribution on 
the wing chenges creating a new starting vortex to 
be dewnstreea. 

This process continues and by the end of tha 
first tlae step, at tn rat, a freo-vortex strip, 
consisting of triangular frae-vertax panels attached 
to the separation edges. Is created. The free- 
vortex strip obeys the conditions given by equa- 
tions (43) and (44) end Its upstreaa edge satisfies 
Kutta condition at the separation edges, equation . 
(45). 

b. At t"tgrat, one needs 7(7, t rat) to deteralne 
the width of the vortex strip. "Since this velocity 
Is unknown so far, a good estlaate U taken as the 
velocity at the preceding tlae step 7(r,tg), which 
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It coaplttily knoMi frt» tta solution corrotpendlnt 
to tho Initial condition. Hlthln tho tlao tup At. 
tho dltpUctaont aF of tnjr ponol nodo It found froa 

A? • ♦ At) - 7(t^) • At • TCr.tjj) (47) 


7(r,t^) • -0, ♦ t»(7,t^) 

(«*o for tiM contidtrod coto) (4t) 

^Cr.to) • (49) 

k-1 


7|^(Fit^) It pivoti Oy oouotlon (19) (SO) 

and It too total mabor of ponolt on tho winp. 


At t>t Ht, OM alto noodt tto uortlelty of 
tho trlangwTarjfortox panolt foralno tho ^roo- 
verut strip, u^Cr.t.nt). POr oach triangular 
panol, M axpritt Its fivo unknown coofficlintt 
(a,)-(a ). datcrihing Its llnoar vorticity dittrl- 
hut1on,*1n tom of tho fivo unknown coofficlontt 
of tho adjacont beund>vortox panol, at tho topora- 
tion odgot, at t-tg (a.-aK of tho bound panol art 
atroady known) and at (n’^At (a, -a. of tho bound 
panol aro still unknown). This it achlovad by 
Mtitfying tho following conditions. At tat.4«t, 
tho vorticity It. contimious at tho global nmt (on 
tho toporatlon odgo) botwoon tho bound- and f^o- 
vortox panolt. At t«to^At. tho fluid porticlos 
along tho downttroaa odgot of tbo triangular panolt 
aro tho taao particlot which oxittad at tbo odgo of 
tmratlon at tng. According to Holholtz thooroo 

■ G»*f)v]. tho vorticity of thoto porticlos 


changot as thoy aro convoctod downttroaa. lut 
Kcrdl., to Ut.l. thvir. elrc.1.. 

tion around thoso particlot rawlnt constant and 
honca additional orations art written to utitfy 
Kalvin thooroo botwoon t^t^ and t-t^^At. 


_ Noxt. tho vorticity distribution on tho wing 
w(r,tgtAt) It obtained froa tho overdotaralnod sot 
of algebraic equations. 


c. At t*tg>2At, a new froo-vortax strip It created 
along the separation odgot. tho first shod free- 
vortex strip It convoctod downttroaa under the 
condition ^ , g* <tep (b) It repeated to find 

^he locatiwis of tho freo-panels nodes and 
»(r,tg^2At). 

d. The steady stau It reached once tho change In 
tho vorticity distribution It lost than a pro- 
tcrlbod error. 


A.Cp(r,tk) • - JYj(r,t^) • C’<rr,t,j) ♦ (o. - lExr)] 
»A»(r,t,^) 


It 


( 91 ) 


where 7.(r,t^) It the Jiau In the ungentlal velo- 
city acutt the wing, which It coapletely known 
frTgt the vorticity vector at this loutlon, _ 
«(r,t^l. Tho first tora In the brackets bt(r,t^) 
It uTculatad froa equation (19). The last tora 
on tho right-hand tide It ulculatod at follows: 


»a4(r,t^) Ad(r,t^) - AaCf.t,^^,) 

si • at ^ 


r(r,t^) - r(f,t|j,,) 

At 


(92) 


Y. CoMutor frodTia 

A coaputar pregraa It dovelopod to liuldaant 
tho aothods of solution of tho ttaady and unttaady 
flow probloas. Tho pregraa It divided Into throe 
aajer parts and each part coMittt of tovoral tub- 
routlnot. 

The first part dealt with tho wing gooaotry, 
wing panelling (bound-vertex panolt) and woke 
panelling (froo-vortax panolt). Uuoa of Initial 
wake panelling dopondt on tho Initial conditions 
of tho unsteady arablaa. No wake panelling it con- 
sidered If tho unsteady problaa tUrtt froa rest. 
For tho unttaady problaa with staady-ttaU Initial 
condition, wake panelling It considered. Sovorol 
tubniotlnot aro used to gonorau tho global and 
local coordinates of tho panolt and tho direction 
cetinot of tho local ‘aerdInaU syttoat. 

Tho second part dealt with tho Initial condi- 
tions. For the unsteady problaa with tteady-ttaU 
Initial condition, tho steady-flow problaa It 
Initially solved. Tho boundary conditions on the 
bound- and froo-vortax panels are utltflod (Sec- 
£lon III. 2. 3.) through several tubroutinos. For 
«(P), thoto tubroutinos fill a rectangular aatrlx 
[N.. (nuabor of boundary-condition equations) x 
' N^*tn«iBbor of undotoralnod coofficlontt of tho 
bound-vorticity distribution)] with the coeffi- 
cients of tho boundary-condition aquations; no- 
ponetratlon conditions, continuity of vorticity 
conditions, Kutta conditions, and tyaaotry condi- 
tions. T)m overdotaralnod tot of equat1ons_,1s then 
solved using a least-square solver. Per U(r), tho 
boundary conditions on tho froo-vortax panels aro 
enforced using a separata subroutine. Once tho 
boundary conditions aro satisfied, tho not-prossure 
coofficlontt are calculated as given In Section 
III. 2.4. If tho Initial condition Is rest, tho 
problaa It solved without any froo-vortax panels. 


IV. 2.2. Pressure Distribution 

To calculate the net-pressure coefficient at 
any point on the wing surface r at any tiae 
U ■ t k . At, we apply the unsteady lornoulirs 
e^tlsn 


The third part of tho prograa daals with tho 
problai at t > tg. Several subroutines are used 
to satisfy the boundary conditions (Section 
IV. 2.1.) at t;ta subsequent tine steps t|(-tgH*t; 
k>1,2,... One subroutine is used to generate a 
new shed free-vortax strip and another subroutine 
Is used to convect the previously shod free-vortax 
strips. The vorticity distribution of tho panels 
forolng the previously shod free-vortex strips 


f 
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^(r,t|,) Is found froi inothtr subroutliM wdleli 
•nforcos Ksivin thooroa. Tht unknown vortlclty 
coofficloflU of tho ponols fonolnf tiM nowly shod 
vortox«slMot Sir Id of* rolatod to tlioso cooffi* 
clonts of tito wing ptnoli oxlstlng ot tlw sopor • 
tion odgos. 

Othor suOroutInos or« ustd to fill tho lattrlo 
rtprtsontlng tlto wing Ooundory conditions^ no> 
ponotratlen condition, continuity of vortlclty 
condition, unstoody Kutto condition and syaaotry 
condition, TIm loast-soMro solvor Is usod to 
obtain ^ vortlclty distribution at tbis tiM 
stop, Z(r»x^^). 

Tbo unstaady lomoulirs oduatlM, tguatlen 
(II), Is tbon usod tm^gh a saparato subroutino 
to calculato tbo not*pi^suro eoofficlont. 

Flguro 7 sboMB a flow cbort of tbo coaputor pro- 
graa. 

11. Wuatrlcal Enaaolos 


11.}. Stoidy now 

Tfw dovolopod c oap u tor prograa Is usod to solvo 
for tho staady flow past a roctangular wing having 
sldo-odgo soMratlon. Tbo wing Is of aspoct ratio 
of ont at I./* angli of attack. Tbo wing Is 
dividod Into Ixl guadrllatoral bound-vortai panols, 
tbo trailing adgo froo-vortojt shoot Is dividod 
Into Ixl guadriiatoral froa-vortax panols, and 
tho sido adgo froa- v ortoa shoot Is dividod Into 
I vortox strlpB and oach Is dividod Into a dif- 
foront mabor of guadrllatoral froo-vortox panols 
such that tho last panol In oach strip occupy tbo 
saao cbordwiso station as that of tbo last panol 
of tbo trailing adgo f r oo- v o r tox panol. 

figuro I shows a typical co n vargad solution of 
tbo spanatlso and cordwlso coapononts of vortlclty 
at tho local nodos of tho guadrllatoral vortox 
panols. At any nodo, tho upper mabor Is tho span- 
wiso conponont w, and tho lowor mabor Is tho cord- 
wlso coi^nont «£. It Is soon that tho continuity 
of vortlclty condition, Kutta condition and tho 
syaaotry condition art utisflod at tho coHon 
nodos. at tho nodos of tho trailing and sido odgos, 
and at tho notas of che lino of syaaotry*, rospoc- 
tivoly. 

figuro 9 shews tho spanwiso variation of tho 
section noroal-forca eoofficlont at throe 
Itoratlon stops. Tbo solution canvorgos In tho 
third Itofstlon stop and Is In good aa r oasio n t with 
tbo axporloontal dau of Seholtx'*. tho figuro 
also shows tho solution of tbo saao caso obtalnod 
by tbo NOV oo tb o d with Ixi bound-vortox linos. 

Ono can concludo that tbo WW-awtbed undorostl- 
■stos tbo noroal-forca eoofficlont near tho wing 
tip. If tbo mabor of botad-vgrtox linos of tbo 
NOV o ot hod Is Incroasod to 9x7*, tho solution 
agroos with that obtalnod hy tho NNVHsathod with 
Ixl bound-vortox panols. This cloarly shews that 
a loss nurtor of vortox panols olvos tho saaa 
accuracy at that obtalnod hy a largo nuobor of 
bound-vortox linos. 

Figuro 10 shows tho convorgod solution of tho 
tponwlio variation of tho suction pitching-oeaont 
eoofficlont for tho wing. Tho rosults of tho WIV- 
oothod with ftxi and 9x7 bound-vortox linos aro 


also Includad In tho figuro. 

Figuro 11 shows tho spanwiso variation of tho 
not-^ssuro eoofficlont at difforont chord stations 
with difforont nuaibor of panols for a roctangular 
wing at IS* anglo of attack, Tho eorrospending 
rosults of tho NOVHBothed aro also shewn In tho 
figuro. 

Tho prosont eooputar iN'ograai :« eonputatlonally 
officiant when It Is coops rod with othor existing 
codas which use higb-ordor doublot distribution. 

Tbo CPU tiM on tho CVICb 171 tbit caso Is about 
200 seconds. 

V I .2. Unstoody Flow 

Next, the dovolopod coopwtor pmorao Is used to 
solve for tbo unstoody flow past lopulsivoly started 
wings froai rest without side odgo separation. In 
tbs two axasmlos, tbo wing Is dIvIM Into IxS 
guadrllatorai bound-vortox panels with a sine dis- 
tribution In tbo cbordwiso direction and a cosine 
distribution In tbo spomrito direction. In tho 
proswt COSOS, tho dlsMoslenloss tlaw stop Is 
oqui volant to 0.49 whilo tbo root-chord length Is 
5 units. 

Fimro 12 shews tho distribution of tho lift 
eeoffleiMt for a roctangular wing of aspect ratio 
of tluroo at 5* angle of attarfc for t«2 and t«22. 

Tho prosont rMults nro coaparod with tho etoody 
iMBorIcal data of rtfertneo (19) w h o re 1M panols 
of constant potantlal function aro usod. It Is 
also coaparod with tbo oxpor1a«ita1 data of 
Roferonce (20). Altbou^ wo used 21 panels in tbo 
prosont ease, tbo rosults coapsros wall with tbo 
othor Huaorlcal and oxporlaontal data. 

Figuro 13 coaparos tho growth of Indlelal lift 
for tbo saao wli^ cans1d«*od above with tho mosorl- 
ca1 data of roforonco (17) whoro 100 panels of 
constant doublets aro used. 

Although wo usod 29 panols In tho prosont caso 
tbo rosults coaporo wall with tho ether nuaorical 
data. 

m 

Figure 14 shews tho dlstr'i/utioa of not pres- 
sure eoofficlont for a roctsi»iular wing of aspoct 
ratio of tae at 20.9* angle of attack for t>2. 

10, 21 and 22. On tho CV4CK 179 coaputor, tho CPU 
tiao fOr each caaa with 9x9 bound panels and 22 
tloa stops Is about 10 ninutas. 

Currently, work Is underway to Incroaso tho 
moibor of panols, uso tho officiant far field 
calculations, and calculate cases with side- and 
loading-odgo separations. 

VII. Concludlno Nworfcs 

Steady and unatoody Nonlinear Hybrid Vertex 
Nothed for lew aspoct ratio wings at largo angles 
of attack Is dovolopod. Tho aothod uses vortox 
panels with first-order vortlclty distribution 
(ogulvalont to soeond-erdor doublet distribution) 

In tho near field calculations. In tho far-flold, 
tho distribution vortlclty Is rotecod to concen- 
tratad vortox lines whoro tho slnplor llot-Savart's 
law is ai^leyod for tho volecity-flold calculations. 
Tho oothod Is applied to steady and unsteady flow 
preblans without any restrictions on tho order of 
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Mgtiltua* of the dlsturboncos In th« flow field. 
The presMted mjaerlcel results shew thet the 
eethod pro<hices eccurate results and It Is coa- 
putatlonally efficient. 
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-F No-P«n«trit1on Condition. Eq. (29) 

O Continuity of Vorticity Condition, Eq. (30) 

• Kutta Condition, Eq. (33) 

Q Sj«Mtry Condition, Eq. (34) 

* Kinontic and Oynaalc Conditions, Eq. (35) 

Figura 5. Dttalls of tha Boundary Conditions fOr 
tba Staady-Flow Problaa. 




Figura 6. Typical Sourcts of Unsteadiness for tha Figure 7. Flow Chart of the Computer Program For 
Unsteady Symmetric Flow Problem tha Steady and Unsteady (Tima Domain 

Approach) Flow Problems. 
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Spanwisa Stations for a 
Ractangulir Ming, AR*3, a*5* 

SxS Bound Fanals, zaro thicknassi 
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Figura 13. firowth of Indiclal Lift for a 
Ractangular Wing, Aft«3, a*5*, 
5x5 Bound Panals, zaro Thlck- 
nast. No Tip Saparatlon. 



prassura coaffldant at 
diffarant chordal sa station 
for a Ractangular Wing, 
AR-2. a*20.5*, 5x5 Bound 
Fanals, Zaro-Thicknass, No 
Tip Saparatlon. 
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RECENT IMraOVEMENTS IN THE PREDICTION OF THE LEADING AND TRAILING 
EDGE VORTEX CORES OF DELTA MINGS 

OtMW A. Kandll* «nd Lakslmnan Balakrlthnan** 
Oapartaitnt of Nachanleal Englnaarlng and Ntehanics 
Old Dominion Univarsity, Norfolk, VA 


Abstract 

Tht rKantly Hodiflad Nonllnaar 31scrata« 
Vortax (MNDV) SMthod has shown a ramarkabla SMCcass 
In pradicting, for tha first tloa, tha latast ax* 
parlmantal data publlshad by Hummal on vortax for* 
matign bahind a slandar dalta wing at an ang1a of 
20 . 5 ®. 


turbulant boundary*! ayar flow, such a paak Is 
hardly noticad and tha prassura paak corrasponding 
to tha primary-vortax cora Is largar than that of 
tha laminar boundary*! ayar f!ow*. A third typa of 
flow Involving a "toi-tlary" saparatlon may occur 
batwaan tha llnas of primary and sKondary sapa* 
rations dua to tha advarsa spanwlsa prassura gra* 
diant ganaratad by tha sacondary*vortax cora^iS. 


This papar prasants tha racant davalopmants 
In tha NNOV«mathod in ordar to accurataly pradlct 
tha location of tha trail Ing-adga vortax cora and 
tha surfaca prassura distribution. 

Noraovar, a Viscous Cora modal, basad on tha 
boundary-1 ayar* ilka approximations (quas1*cy11n* 
drical approximations). Is prasantad. Tha ra- 
sultlng parabolic aquations, with tha outar-adga 
boundary conditions obtalnad from tha Inviscid* 
modal solution using tha NNOV*mathcd, ara Inta* 
gratad using a finita-diffaranca marching tach* 
nlqua. Typical valocity profllas of tha loading 
and trailing vortax coras hava baan usad to 
Initlata tha stap*by*stap marching tachniqua. 


I . Background 

1.1 Existing Exoarlmantal Data and Concluding 
Ramarks 

Tha saparatad fraa*shaar layars amanating 
from tha laading adgas of highly swaptback wings 
roll up spirally Into two countar rotating vortax 
coras. Tha vortax coras ara fad, through tha 
fraa-shaar layars, with vorticity shad along tha 
wing loading adgas from tha boundary*! ayar flow 
on tha wing surfacas. This flow saparatlon Is 
known as tha 'primary saparatlon* and It has a 
dominant offset on tha wing aaredynamlc charac- 
tarlstlcs duo to tha largo strength of Its vortax 
coras and thair proximity to tha upper surface 
under tha primary vortax coru. 

Tha advarsa spanwlsa pressure gradient dua 
to tha primary vortax cora Induces ’secondary 
separation* of tha boundary-layer flow on the 
upper surfaca. Oapanding on the angle of attack 
and tha wing aspect ratio, tha secondary sapa* 
rated flow forms either an additional fraa-shaar 
layer or a bubble. In the range of moderate to 
large angles of attack, a secondary freo*sh«ar 
layer davtiops, rolls up spirally in an opposite 
sansa to that of tha primary fraa-shaar, and forms 
a sacondary-vortax cora of much small at strength 
than that of tha primary-vortax core 

Oapanding on the typa of boundary-layer flow 
on the upper surfaca, tha sacon^ry-vortex cora 
may affact tha prassura dlstrlFutlon on this sur* 
faca. For a laminar boundary-layer flow, a small 
prassura peak Is produced batwaan tha secondary 
and primary Unas of separation while for a 

^ssoclata Professor, Mambar AIAA 
**Graduate Student, Hember AIAA 


It Is seen from tha description given above 
that secondary and tertiary separations ara viscous 
phenomena and cannot be modeled by using Inviscid 
analysis only. Howavar, thair affects ara small 
particularly whan turbulent boundary* layer flow 
exists on tha upper surfaca of tha wing and hanea 
they are neglected. Tharafora. whan Inviscid 
analysis Is usad to pradlct prassura distributions, 
ona tmjst compart the surface pressure distributions 
with those axparlmantal data corresponding to 
turbulent boundary-layer flows®. 

The fraa*shaar layer amanating from tha 
trailing adga Is of opposite strength to that of 
tha primary fraa-shaar layer. Within a short 
distance bahind tha trailing adga (a distance of 
a 1/4 root-chord lanoth for a dalta wing of aspect 
ratio of one at 20.5® angle of attack®), tha trail- 
ing fraa-shaar layer rolls up spirally In an oppo- 
site sansa to that of tha primary fraa-shaar layer 
and forms two countar rotating vortax coras. Tha 
trail Ing-adga vortax cora has tha same sansa of 
rotation and almost tha samO spanwlsa location as 
those of tha secondary vortax cora, although each 
originates from a different phanomanon. The formar 
Is due to an adga saparatlon while tba latter Is 
dua to a surfaca saparatlon. Hummel® concluded 
from his latast axparlmantal maasurments for tur- 
bulent boundary layars at tha trailing adga that 
the sacondary-vortax cora decays rapidly bahind tha 
trailing edge and Its remains roll up Into tha 
trailing vortax. 

Figure 1 (rap.'oducad from Rafaranca 5) shows 
tha vertex formation bahind a slender delta wing. 
Figure lb corrects tha location of tha trail Ing- 
adga core as given earlier'** by Figure la. It 
also shows that tha sacondary-vortax cora has the 
same spanwlsa location as that of tha trail Ing-adga 
vortax cora but It does not show thair Interaction 
which Is still unknown. 

Predicting tha correct formation of tha trail - 
Ing-aoge vortex cora Is assantlal for predicting 
the correct aerodynamic loads near tha wing trail- 
ing edge due to Its proximity to this edge. More- 
over, a correct prediction Is of paramount Impor- 
tance to probimns Involving high-angla-of-attack 
aerodynamics for tha canard-wing configurations. 
Consequently, ona concludes that any prediction 
method which employs a prescribed shape of the 
trail Ing-adga vortex sheet (the so-called 'fixed 
design wake'l®* will be strictly limited to 
Isolated-wing problems away from tha tralling-edge 
region. 

Tha size of tha primary-vortax cora and 
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vliceilty of tht flow within « narrow rtglon around 
tha cert ctnttrlint play luportant rolti In tht 
atrodyntalc chtrtcttrlttlca of low-asptct ratio 
wings at largo angits of attack, for thatt typtt 
of wings • it Is wall known that tha lift coaffi* 
clant Incraasas nonllnaarly with tht Incraasa In 
tha angla of attack. This Is attributed to the 
Increase of the strength of tht prlsiary vortex 
core. However, the Increase In tha lift cooffi* 
clant Is llnltad by tN occurrtnce of vortax>cora 
breakdown ovar tha wing and henca there exists a 
maxlMai angle of attack corresponding to the smxN 
lauai lift. 

The phenoMnen of vortex -core breakdown Is 
defined as an abrupt Increase In tha cere diaaieter. 
It occurs due to the adverse sw1r1;1nduced pressure 
gradients along tha core ax1s'^» Using the 
radial coaponent of tha aoaientua aguatlon of tha 
boundary>1ayar Ilka aguatlons (for guat1«cy11n* 
drical vortex cores), one can show that an Increase 
In the swirl velocity of the flow, due to an 1n« 
crease In the core strength. Increases tha axial 
pressure gradient along tha core axis. Therefore, 
tha flow near the core axis loses Its axial aoeMn* 
tua and swelling of the vortex core develops. 

For aoderate angles of attack, vortex «core 
breakdown develops far downstreaa froa the wing 
trailing edge. As the angle of attack Increases, 
vortax«core strength Increases and the breakdown 
point travels upstreaa toward the trailing edge. 

As the breakdown point approaches the trailing 
edge, the slope of the lift versus angle of attack 
curve decreases until the siaxlaua lift Is reached. 
Thereafter, the lift coefficient drops sharply, 
the breakdown point crosses the trailing edge and 
the wing stalls. 

To account for all these changes, one has to 
Include the viscosity within the vortex core In 
order to have a realistic model In this region. 


I.2.Ex1st1ne Invlscld Hodels and Concluding Remarks 

In most of the existing Invlscld analyses, 
primery flow separation Is only considered. Tha 
attached boundary layers are replaced with bound- 
vortex sheets while the free-shear layers are re- 
placed with free-vortex sheets. Accordingly, the 
Invlscld region representing the vortex core should 
be represented by a conical spiral vortex sheet 
which ends with a concentrated vortex line along 
the core centerline. Although this theoretical 
representation of the vortex-core Is Ideal for the 
Invlscld eedel, the usage of many turns within the 
core region creates nunerlcal problems In a three- 
dimensional model which does not assume slender- 
body approximations. The numerical problems arise 
due to the proximity of the vortex layers of the 
spiral and due to the large number of small -size 
panels needed to accurately model the turns. How- 
ever. It was found, according to the numerical 
results and their excellent comparison with the 
experimental data, that only one turn of the spiral 
need to be accurately modeled while the remaining 
turns of the spiral are replaced by a cut ending 
with a concentrated vortex line along Its edge. 

In all the Invlscld models, the wing edges, 
where separation occurs, are assumed to be sherp 
so that the separation lines are known a priori 
and hence viscous modeling Is not needed to deter- 


mine the lines of stparatlons. However, for wings 
with round edges, viscous modeling Is needed to 
determine these lines'*. Moreover, vortex-break- 
down points are assumed to be far downstream to 
that tha variations In the size of the primary- 
vortex core can be neglected. In fact, this as- 
sumption limits the large angles of attack at 
which Invlscld modeling Is applicable. With the 
vorticity confined to the Inner boundary of the 
flow region (bound- and fret-vortex sheets), the 
resulting model Is a potential -flow model governed 
by Laplace's eguatlon end certain boundary con- 
ditions'*. 

The existing Invlscld models can be divided 
Into four main groups. Next, each group Is pre- 
sented and evaluated. 

Tha first group of models uses slender body 
and conical flow assumptions'** Modeling of 
the primary-vortex core and Its feadipg sheet was 
first Introduced by Irown and Michael '. However, 
the feeding sheet In their modal was taken as a 
planar surface and henca It did not represent the 
rnal rollad-up vertex sheet. Mangier and Smith'* 
Introduced the first realistic model of the primary- 
vortex cere and its feeding sheet. However, this 
model does not account for the feeding-sheet defor- 
mation In the chordwise direction duo to the 
slander-body assumption. These models satisfac- 
torily predict the pressure distribution on the 
front portion of the wing surface. In the rear 
portion, the models fall to predict satisfactory 
pressure distributions because Kutta condition can- 
not be satisfied at the trail Ina edge. Such models 
were limited to slender delta pianforms. 

The second group of models uses the old NOV- 
method<*'3h . The most obvious drawback of the old 
model Is the lack of a realistic model of the 
primary-vortex core and Its feeding sheet. The 
primary-vortex sheet was modeled by a system of 
segmented vortex lines which were rendered force- 
free during the course of solution and no attempt 
was made to lump these lines Into a concentrated 
primary-vortex core. Although the agreement 
between the calculated and experimental total aero- 
dynamic loads was excellent, the agreement between 
the calculated and experimental pressure coeffi- 
cients was lass than satisfactory In particular for 
highly sweptback wings. 

The third group of models employs doublet- 
panels'^' "• in this method, the wing and 

Its free-vortex sheets are divided Into networks of 
quadrilateral panels. Each panel of the networks 
representing the wing has a biquadratic local 
doublet distribution and a blllnenr local source 
distribution. The panels of networks representing 
the free-vortex sheets have biquadratic local 
doublet distributions. Source and doublet splines 
are used to express the distributions of singu- 
larities on the networks In terms of discrete velues 
of singularity strength at certain standard points 
on each network. The boundary conditions and con- 
tinuity of singularity strengths across abutting 
networks are enforced at certain standard Mlnp on 
each network. The results of this method*** *' are 
generally good when the solution converges. Ap- 
parently, the difficulty In obtaining convergence 
Is due to the failure In satisfying the continuity 
of the derivatives of the doublet strength across 
abutting networks. This Is equivalent to the exis- 
tence of concentrated vortex lines between abutting 
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nttworki. 

42 

Tht doubitt ptntl method wit utindtd to 
Includo tht ifftct of intrilnmnt of tht prlnwr/ 
vorttx eertt through in mptrlcil ipproich. Tht 
rtiuUt Indlettod thit tht tntrilnmtnt Incrtiitd 
tht normiUferct coefficient tubttintlilly over 
the txptrlmtntil vilutt. 

Thli group of model t do not iccount for the 
chordwlit deformtlon of the pr1mery<verttx sheet. 
Centeguently, the prlmery vortex sheet einnot be 
fed three*d1meni1oni11y with the shed vorticity. 
Moreover, the trill 1ng>edge free>vortex sheet wet 
represented by i "fixed design wtke." 

Therefore, none of the models given In the 
first ind third groups Is cipible of predicting 
the deformtlon ind Interictlon of the trill Ing- 
edge vortex sheet ind Its vortex core. Addltlonil 
references of the three groups given ibeve cm be 
found In references 43 ind 44. 

The fourth group of models Mploys i nonllneir 
hybrid vortex method'^t w, in this method, con* 
tlnuous-vortlclty tnd vortex*11ne representitlons 
of the wing ind Its sepirited free-sheir liyers ' 
ire used. Continuous vorticity Is used In the 
neir>f1e1d cilculitlons while discrete vortex*11nes 
ire used In the fir*f1e1d cilculitlons. 

The wing ind Its free*sheir liyers ire divided 
Into plinir quidrlliteril pinels ntving f1rst*order 
vorticity distribution. The lerodynimlc boundery 
conditions ind continuity of the vorticity dlstrl* 
butlons ire Imposed it certiln nodil points on the 
pinels. To sitlsfy these conditions, in Itentlve 
technique Is folVowed which iltemitlvely yields 
the locil vorticity distribution on the bound* 
vortex pinels end the shipe of the free-vortex 
pinels. This method his been used to cilculite 
the steidy distributed end totil loids on plimr* 
1ow*ispect*rit1o rectengulir wings. The results 
hive shown thit the spinwlse virlitlons of the 
loid coefficients ire In good igreement with the 
experlmentil d&ti. Compirlsons of the results 
with those of the N0V*method hive shown thit the 
hybrid method requires less number of vortex pinels 
for the sime iccuncy. Currently the method Is 
extended to Include lending ind side-edge sepi* 
ritlons for steidy*flow problems end for unsteidy- 
flow problems due to smill osclllitlons of wings 
iround lirge mein ingles of ittick (flutter ippll* 
citlons) tnd due to generil nonstitloniry wing 
motions (dynemlcil ippllcitlons). Uork Is underwty 
to use this technlinie to predict the detills of the 
flow field Including the prlmiry tnd trill Ing-edge 
vortex core Interictlons. 

Although the old NOV-method. used In the 
second group of models, wis ilmost tbtndoned In 
predicting the flow dettHs tnd the distributed 
terodyntmic chiricterlstlcs, the modified NOV- 
method pinpointed tnd cured mest.ef the problmn 
encountered with the old mthod^. It enjoyed t 
remirkible success In predicting, for the first, 
time, the litest experlmentil results of Hummel^. 


1.3. Existing Viscous-Core Hodels ind Concluding 
RiBiirks' 

In this piper, we only consider existing vis- 
cous-core models which idopt the quis1-cy11ndr1ci1 


ipproxlmetlons. In the quisl-cyllndrlcil tpproxl- 
Mtlons, we tssurne thit the ixlil gndlents ire 
smell compered with the rtdiil gndlents tnd thit 
tht streimsufices be cylindrictl. These boundiry- 
Ityer-llke issumptlons simplify Nivler-Stokes equi- 
tlons end reduce then to t ptrtbollc set of equi- 
tlons. Approxlnitlons of.thls type hive bgen used 
by 6irtjhore47, hi 1I*«.49, lottel'^. Rut*® tnd 
Milton’' to predict the locitlons of the vertex- 
core breekdown. Girts here tnd Hill used the dif- 
ferential form of the equitlent while Bessel, Rut 
end Hilton used the Int^rel fern ef the equitlent. 

The nusMrIetl results ebtilned from the quisl- 
cylindrieil equitlent hive shewn good igreement 
with the experlmentil ditt pirtlculirly for the 
vortex core behivlor. for the locitlons of the 
bretkdown points, the difference between the pre- 
dicted end meisured locitlons wis within two dli- 
meters of the viscous core. 

In this piper, we use the quis1-cyl1ndr1ct1 
ipproxlmitlont for the nirrow viscous region iround 
the core centerline. Outside the viscous region, 
one can treat the flow as a potentleT flow only or 
as a rotational Inviscid flow followed by in outer 
potential flow. Hence two models with various 
degrees of accuracy emerge. In the first model , 
the vortex core Is represented ty a nirrow ix1 sym- 
metric viscous region iround the core axis. This 
viscous core Is fed with vorticity shed from the 
wing edges through a vortex sheet with only one 
turn and a cut. The flow outside the viscous core 
and the vorcex sheet Is a potential flow. Figure 
2.1. shows 1 tchwitle of this model. In the second 
model, the vortex core region Is divided Into an 
Inner narrow viscous part and in outer Inviscid and 
rotational part. In each part the r^ow Is ixlsym- 
metrlc. Outside the core region, a potential flow 
Is assumed. Figure 2.b. shows a schmtic of this 
model, n this paper, only the first model 1.s 
considered. 


II. Steady Inviscid Model ind Prediction of the 
Inv^sc^d Cores 


II. 1. 

Nomenclature 


AR 

wing aspect ratio 


b 

wing half span 


b(x) 

local half span 


'p 

static pressure 



net surface pressure 


'r 

wing root chord 


ITER 

Iteration Number 


LE 

Leading Edge 


LEC 

Leading Edge core 


TE 

Trailing Edge 


TEC 

Trailing Edge core 


xyi 

wing-fixed coerdlnates, origin at wing 

iyi 

vertex, x-ixis Is along thi root 
y-ax1s Is perpendicular to wing 

chord. 

wind-fixed coordinates, origin at 

the 


trailing edge, x-ax1s Is parallit 

to the 

0. 

freestreim direction 
freestream velocity 


V 

velocity at any field point 



.Vscomponents of V In the wind coordinate 
system 


component of V In planes x • constant, 
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♦ V?)** 

a «ng1t of ttUck 

dlMfulonUts coordlnitit, 
5«I/b, n*y/6» ;*i/b 


U.Z. Culttina N0».wthod_ ind Iti Drwbickt 

In this itctlon. « critical evaluation of tha 
txlatlng N0V*Nethod It given In order to pinpoint 
ita draebackt. In this nethod, the bound •vortex 
sheet It replaced by a bound-vortex lattice and 
the free-vortex iheeti are replaced by lagtMnted 
free- vortex lines. 

Starting with an Initial west for the In- 
clination of the free-vortex tines, the no-pene- 
tratlon condition Is satisfied at certain control 
points on the bound-vortex lattice. The resulting 
set of linear algebraic equations Is solved for 
the circulation distribution. Next, with the 
known circulation distribution, the positions of 
the vortex segMnts of the free-vortex lines are 
calculated by slaultaneously enforcing the no- 
penetration and no-pressure Jump conditions at 
certain control points on these segaients. These 
two steps of calculations represent one Iterative 
cycle. Several Iterative cycles are performed 
until the positions of the free-vortex segments 
or the circulation distributions converge. 

Figure 3 shows a typical converged solution 
of the system of free-vortex lines In three views 
for a delta wing of aspect ratio of unity and IS 
angle of attack. The plan view also shows the 
arrangement of bound- vortex lattice. In the three 
dimensional view the loading-edge core (LEC) Is 
shown. This core Is calculated aftor the solution 
converges and it represents the centroid of the 
leading-edge vortex system. 

Comparison of height and spanwlse position of 
the calculated centrald with those of tha measured 
leading-edge core^ Is given In figure 4. Al- 
though this comparison was encouraging, the calcu- • 
lated centroid does not model the physical vortex 
core where tha vortex core Is continuously fed 
with forticity from the leading edge through the 
free-sheer layers. Using the system of free- 
vortex lines, the total -aerodynamic loads were 
accurately calculated but the calculated surface 
pressure distribution wet unsatisfactory. On tha 
other hand, using the centroid of the leading-edge 
vortex system instead of the leading-edge vortex- 
system Itfolf. the calculated surface pressure 
Improved^B (this does not satisfy the no-penetra- 
tion condition on the wing since the centroid Is 
calculated after the solution converges.) It 
should also be noted here that replacing the 
luding-edge vortex system by Its centnid Is 
similar to tha model used by Icgandre'^. 

Figure 5 shows a recent converged solution of 
tha system of free-vortex lines with a long de- 
formed wake. It can be teen that tha trace of the 
tra 1 1 1 ng-edge vortex system In cross planet Indi- 
cates that the sheet tends to deform upwards 
showing a tendency to form a trail Ing-edge vortex 
core. However, the cross-flow planet taken further 
downstream show that the free-vortex lines leap- 
frog. This does not represent the real flow. 

Figure 6 shows four of this cross-flow planes 
taken perpendicular to the wind direction. 


It It clearly teen from the few examples 
given above, that the existing model of the NOV- 
method does net realistically medal the leading- 
end trail Ing-edge vortex cores. Therefore, the 
model and the numerical technique mutt be modified 
In order to obtain realistic vortex-core modeling. 


11.3. The WIOV-Hethed and Model 1.w of the Primary- 
Vortex tore 

In tha MMM-Methed a realistic modeling of '"'t 
primary-vortex core region It Introduced. The oid 
NOV-method Is only used In the first Iterative cycle 
to initiate tha roll-up process and to calculate 
the centroid of the leading-edge vortex system. 

Next, the leading- and tra111ng-edge vortex segments 
are replaced by smaller segments. Then, the Itera- 
tive cycles proceed. 

In a typical Iterative cycle, each free-vortex 
line of the leading-edge vortex system Is allowed 
three-dimensional ly to rotate around the most 
recently calculated centroid a prescribed portion 
of a turn (1/4, 1/2, 3/4, or 1 turn). This Is 
done by continuously checking the coordinates of the 
free-vortex segments (as they are adjusted) with the 
location of the calculated centroid. Once this Is 
achieved, the remaining free-vortex sefnents of each 
line are dumped Into the calculated centroldal line 
3f the vortex system. The Iterative technique Is 
followed until the circulation distribution con- 
verges . 

Figure 7 shows typical solutions at different 
Iterative cycles. It can be teen that after one 
Iterative cycle, the system of free-vortex lines 
shows good roll-up. The converged solution. Indi- 
cated by ITERwfi, shows the leading-edge core and 
Its feeding free-vortex lines. It can also be teen 
that the free-vortex lines continues to feed the 
LEC beyond the trailing edge. Tha trail Ing-edge 
core Is also Indicated on the figure. This will 
be clearly seen In the cross-flow planet discussed 
In the next section. 

Figures 8 and 9 show converged solutions >*<ig 
a 12 X 12 bound-vortex lattice for 1/4 and 1/2 turns 
of tha free-vortex lines. 


II. 4. Ntmerlcal Cxynolw and Comparisons with 
HueweVs Experimental feetuTEs 

The calculated circulations around the free- 
vortex lines emanating from the wing trailing edge 
was found to be of opposite sign to that of the 
free-vortex lines emanating from the wing leading 
edge. This Is In agreement with Hiamael's measure- 
ments. The difference In the signs can be explained 
at follows: 

The primary-vortex core creates large suction 
pressure peak and hence the pressure continuously 
rises In the spanwlse direction from the wing axis 
to the location of the suction peak (this Is com- 
pletely opposite to the sn^nwlse variation of 
pressure for wings with 1urge aspect ratio and small 
angles of attack.) Consequently, the circulation of 
the spanwlse bourd-vortex $« 9 nents Increases in the 
spanwlse direction. In order to satisfy the spatial 
conservation of circulation at the nodes of the 
bound-vortex lattice, the spanwlse Increase of 
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circulation rtquint tho circulation of th« chord* 
Witt bound*vorttx it^Ronti to bt of oppotitt tipn. 
In tht prtttnt nodtl, tht 1t«d1ng*tdgt vorttx 
tyttM originatti from ipanw1tt*vorttx itgaitnti at 
tht Ttading tdgt wh11t tht tra111ng*tdgt vorttx 
tyttaai orlglnattt froai chordwitt- vorttx ttgatntt. 

According to tht difftrtnct In sign of tht 
circulation of tht tralling-tdgt vorttx t/ttm froa 
that of tht 1tad1ng*tdgt vorttx lyitt*, ont axptcti 
tht tralling-tdgt vorttx lyatta to roll up In an 
oppotitt ttnta to tht roll-up of tht Itadlng-tdgt 
vorttx lyttaia. 

In rtftrtnca 4€. tht rttultt of tht MNDV- 
tathod Ming 1/4 and 1/2 tuma only for tht Itadlng- 
tdgt vorttx lints wtrt prttanttd. In this ptpar, 
wt add to thttt casts tht rasults of tht 3/4 turn. 
Tht lattst rtfults citarly show tht capability of 
tht MHOW-aitthod la rtproduclnoeand conflrtlng tht 
txptrlawntal results of Huasatl’. 

FlgurtslO-12 show coaiparlsont bttwttn tht 
txptrlaitntal and nusMrlcal rasults for loading- and 
tralling-tdgt vorttx shttts and flow dlrtctlons In 
cross-flow plants ptrptndicular to tht wind dlrtc- 
tlwi. Tht maatrleal results art drawn at tht laaw 
ictlt at that of tht txptr tantal rttultt. Tht 
predicted sizes and locat cis of the leading- and 
tralling-tdgt vorttx shat:' art In axctllant agrtt- 
ntnt with tht Mtsurcl or.^i. Tht rttultt of the 
3/4 turn show tht the correct roll-up of the trail* 
Ing-tdM vorttx sheet and tht correct locations of 
tht LEC and TEC. This Is exactly what wt expected 
when wt Increased the roll-up frot 1/2 to 3/4 turn. 
With tha 3/4 turn of tht Itadlng-tdgt vorttxsysttm, 
the roll-up tlghtans and larger vtlocltlts art 
Induced at the tralling-tdgt vorttx sheet which 
causa It to dtfona upwards and leftwards. Out to 
the large curvature of the trail Ing-edgt vortex 
sheet (Figure 12. 3/4 turn cate), out can see a 
saitll leapfrog at tht TEC location. This Is 
attributed to the Insufficient nutbtr of vorttx 
lints used to todtl this large curvature. 

Figure 13 shows tht results for an aspect 
ratio of 1.43. It Is tttn, by conparlng thttt 
results with the corresponding results for the 
aspect ratio of «ie, that the TEC develops earlier 
than that of the aspect ratio of ont. 

Figure 14 shows cowparlont bttwttn tht pre- 
dicted and measured static-pressure contours In 
different cross-flow plants. The predicted sizes, 
locations, and levels of tht pressure contours are 
In good agreement with those of the measured data. 

Figure 15 shows comparisons of the predicted 
and measured spanwist net surfaca pressure variation 
at difference chordwitt stations. It It teen that 
tht results of tht 3/4 turn are tore accurate than 
that of the 1/2 turn. With 3/4 turn, tht ro11-up 
of the Itadlng-tdgt vorttx systam tightens and 
produces better predicted distribution than that of 
tht 1/2 turn. 

Tht developed compuur program of the MOV- 
Mthod Is computationally tfficitnt. OnaCYBER-175 
■achint, a typical solution using a 12 x 12 lattice 
and Including tht cross-flow plants calculation 
takes 10 minutes o' CFU time. 

Currently the technique Is modified for 
accurate ntar-fltld calculations by replacing the 


concentrated vorttx se^aentt with vorttx panels 
having linear vorticity distributions. Tht vortl city 
functions art expressed In tents of the unknown 
circulations of the original concantrattd vertex 
segments. This Is tht opposite process currently 
used In the Nonlinear Hybrid-Vortex ftathod^St*®. 


111. Steady vitcoM Model of the Vortex Core 
IIM. NomencUture 

A^ grid aspect ratio ■ peMx 

a fraction of the step size In the x-d1rect1on 
b ■ 1 - a 

4c^ characteristic length 

h radial velocity in the computational domain ■ 
w / ; 

k swirl velocity In the compuUtlonal doaialn ■ 
XCv 

m axial station number 

n radial station nuaber 

p pressure 

R ■ b^ 

R Reynold's nuaber . 

r* stretched radial coordinate ■ rR* 

r radial coordinate * 

r (x) stretched outer core radius at any axial 
station 

r,(X|} stretched outer core radius at the Initial 
' axial station 

freestream velocity 

u axial velocity 

V swirl velocity . 

w stretched radial velocity ■ wR* 

ft radial velocity * 

x axial coordinate 

Ax axial step size 

X shape factor ■ r^(x)/r^(*^) 

C radial coordinate In the eomp itatlonal domain 

Eg outer radius In the computational domain 

AE radial step size 


(11.2. Formulation of the Problem 

Starting with tht nondimenslonal Navler-Stokes 
equations for Incompressible axlsymmatrlc steady 
flow, using a boundary-layer-like stretching trans- 
formation and performing an order of magnitude 
analysis, one obuins the nondimenslonal equations 
for quasi -cylindrical vortex flowsA*. 
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This It • parabolic tat of aquatloni anO hanca a 
finita diffaranca narching tacbnigua procaading in 
tha axial dirKtIon la appllcabla to tbit sat. A 
basic diffaranca batwaan tbata aquations and tha 
first ordar boundary* layar aquations It tha 
pratsura tana. In tha boundar>'*1tyar aquations, 
tha pratsura across tha boundary-1ayar thicknatt It 
constant whi la In tha prasant aquations tha pratsura 
It varying In tha radial diractlen at givan by 
aquation (2). 


on only ona boundary; althar on tha axis or on tha 
boundary • const, uhlla tha bout^ry conditions 
on u and Irautt ba spaclflad on tha axis and tha 
outar boundary. 

Along tha cora axis, ; *0, tha folloHing 
boundary conditions ara Impotad: 

. k(x.o) • h(x,0: • 0 (14) 


Tha viscous raglon d1var«t In tha downstraan 
diraetlon dua to tha twirl inducad advarsa pratsura 
graditnt. If tha finita-diffaranca grid it con* 
ttructad In this physical doMln, ona naads to 
fncraata tha ntiabar of grid points In tha radial 
diraetlon at ona procaads downttraa*. For coiapuu* 
tional afficlancy and convanlanca of applying tha 
f1n1ta*d1ffaranca narching tachniqua, It It datir* 
abla to transfom tha physical dlvar^ng donaln 
Into a rectangular conputatlon dosMln. In this way, 
aqual n«aabar of grid points can ba usad In tha 
radial direction at each station of tha axial 
direction (tea figure 16). This Is easily achieved 
by adopting tha trantfomatlon used by Hall^. 

C ■ r(x)/X(x) (5) 


Tha first condition axprastts tha fact that tha 
axial velocity It synnatrlcal about the axis. Tha 
second condition requires tha circulation to vanish 
on tha axlsi l.a. , a rigid body notion exists In 
the Unit at tha axis is approached. Tha third 
condltlwi requires tha radial velocity to vanish In 
tha abtanca of sources or sinks on tha axis. On 
the outar surface, the following boundary conditions 
must ba spaclflad; 

u(x,Cj) • u^(x), k(x.c^) • kg(x), p(x,Cp) • P,(x) 

(15) 

Tha functions u.(x), k.(x), and p (x) ara provided 
by tha 1nv1sc1d”solutl8n of tha oOtar flow obtained 
In Section 11.4. 


X(x) • r^(*)/r,(*^) 


(6) Tha Initial conditions at tha Initial station 

x^ are given by 


Thus, tha outer boundary ( at any axial station 
bacoMS a consUnt 

?o “ • '’0^*1 5 ■ e®"*t*"* 

Noraovar, tha swirl velocity v Is transfonxad Into 
a variable proportional to tha circulation k, while 
tha radial velocity Is scaled by tha shape factor 

k(x,c) ■ nr(x,r) • lev (8) 

h(x,c) • w(x,r)/X (9) 


Using tha transfenaatlon aquations nf aquations 
(5), (6), (8) and (9) Into aquations (1) - (4), wa 
9*t 


>E . 

r? 7? 

X' k^ Jp 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


Equations (10) * (13) are tha governing aquations 
In tha computational dcnaln. Cqijatlons (10) * (U) 
ara first order aquations In h and p; raspactivaly, 
while aquations (12) and (13) ara parabolic 
aquations In u and k; raseactivaiy. Tharafore, the 
bounoary conditions on h and p must ba spaclflad 


u(x^,C) ■ u,(;), k(x^,c) - k,(c) (16) 

Thasa profiles are obtained from axparlmantal data 
or from previous approximate solutions. Aluma* 
tively, they may ba obtained from tha three 
dimensional Viscous flow solution on tha wing 
surface. 


1 1 1. 3. The Finite Diffaranca Tachnloua and 
Quasillnaarltatlon 

In this section, wa develop tha basic dif* 
faranca aquations required to determine a function 
or Its darlvatlvas along tha axial and radial 
dlractlons with sacond>ordar accuracy. Naxt, we 
develop a 'quasi! inaarizatlon tachniqua* to 
llnaarize the nonlinear terms In tha aquations. 


Figure 17 shows a rectangular compuUtlonal 
domain covarad by a grid systam consisting of 
constant x and c Unas. Subscript a rafars to 
tha point number In tha x direction while subscript 
n refers to tha point number In tha c direction. 

Tha developed aquations ara evaluatad at the x 
location which Is defined at follows: ^ 

*c ■ *a '{Vl' *m’ 

where a r b • 1 


Expanding tha functions ^ and 6^^, ^ about point 

(c.n), adding and neglecting terms of 0(ax)^ or 
higher, we get 


-c,n 


®**«,n * ,n 


0(ax)' 


( 18 ) 


Subtracting tha two expansions 
we gat 


*'B,n ***** ®B*I ,n' 


6 



*®c n 
1*^ 




Ax 




1*6 


»x‘ 


♦ 0(4x)^ 


•y 1 ftt 1 n 9 b-«*0(Ax)t tht itcond ttm on tht right 
hACOMt of 0(tx}< ind tho tguotlon can ha truncataa 
thara and wa gat 


ikil • -%ig ♦ o(4x)* 

)x 4X ' 


(19) 


expanding 8^ and about point (c.n), 

auhtractlna and adding tho aguatlona, rj obtain 
tho following aguatlona; raipoctlvaly 


>6 


^ ♦ 0(4«)^ 


(20) 


* mS x O . * ] i^ixD — A 0(0^)^ (21 ) 

>c (4C) 

Substituting ag. (18) Into agi. (20) and (21), wo 
9«t 

*^c.n • V»1 . n*1 ' Sa»1.n»1^ * ^^^.n»1 ' ^.n-1^ 

)C 24( 

♦ C(40^ 

* ^c.n • *^^w»1.nAl ' ^‘Vi.n * *Vl.n-l^ ♦ 

UCT 




28. „ ♦ S_ . ,) 2 

— Sxfl ♦ 0(4C)' 

(4C)^ 


(23) 


Noxt, WO consldar tho "Quaalllnoarlxatlon Tach- 
nlqua.* Equatloni (10) • (12) Includa nonllnoar 
tanas, a.g.; ,au ju ate. which product 

Uj5, uyj-. . . . 

nonllnoar a'gobralc tanas In tha finita diffaronca 
fora. Sines tha solution of nonllnaar actions Is 
t1»a consusilng, ona naads to Hntarlta thasa tarn. 
For this purposa, wo consldar Taylor's axpanslons 
of two functions ^ and ^ about 6^ ^ and 

N ; ravpactlvaly, and fora tha Hnaarlsad axpras* 


111.4. Dlffartya Cauatlons of tha Sorarning 

.. «ln. (c.n.,). 
Using tha diffaronca axprosslens davalopad In tha 
procadtng taction, ona obtains tha following 
aquation; 

'••♦1,n • Vl,n-1 ^ ''o.n.l ’ ^,n^ 

♦ ^(ISjf)([p,^(2„.j).1]«^ ,n - C*v. 

• “si.n.1 ♦ t'’o • ’3 «,4l,n . 

* Vl,nx1 


*r-4i 


(27) 


(28) 


Equation (II) Is avaluatod at point (c,n<t). 
Using tho diffaronca oxprtstlont ^alopod 1<) tho 
procading taction, ona obtains tho following 
aquation: 

**oAl,n-l “ ^OAl.n * *^^o,n " %,n-l^ 

* (2n-3)^ * '‘■,n-l^ * ^'‘•♦l.n * '*»rl.n-l^^ 

(29) 

wharo 

B„ - -1/Ca^ 


l^hgSSSDSJSkill on 

Equation (12) is avaluatad at point (c.n). 
Using tha diffaronca axprrttlons davolopod In tho 
procading taction, ona obtains tho following 
aquation: 


®o*1 ,n **■♦! ,n * Vrl ,n **a,n ^ Si.n **or1 ,n 
^.n «.,n ♦ 


It llnaarliad at follows 

S H ■ a(G . M 'T S M , r 
“c.n c.n ^ ,n a.n ^.n orl.n 

(«*1) 6, - H J » 0(4,)* 

■ •fl Man 

whara i»b/t 

It should ba notad that aos. (18) and (24) 
to dariva aq. (2$). 


» an 

\ ''•♦l.n-l * *n ''•♦l.n * ^a♦l,n-l " ®n 

wharo 

(31) 

(24) 

1 at 


(32) 

It point 

•n ‘ »n - *• «n 

(33) 


Cn • *<«n * "n’ 

(34) 

(29) 

“n * '*f«n * ^n^ 'Ss.ntl * ^"n * **» “n^^.n 


(21) 

■•*^®n * ■»n^‘‘o.n-l 

(35) 

ara utad 



(36) 
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-rV 

- . i ij 


■'« • ‘ Vi,.l 

- ^ * Vl.n* ■ 


whtre 



(47) 


x-Moaentun Equation 

Equation (13) Is tvaluatad at point (c.n). 
Usinq the differtnee expressions developed In the 
preceding section, one obtains the following 
equation: 


^n “n+l,n+l * ®n “lerl.n * T»ri,n-1 ' *^n 


where 

*n • *<®n * ^n^ ^ ®n “w.n 

®n * *1 " ^*®n * ®n^%i,n+l " %,n-l^ 

^ “ *®n - **n - ®n “«.n 

♦ [=*„ * 2M„ - (R*l)(u„,^l - “ai.„li)9n] 


(39) 

(40) 

(41) 




“m.na-1 ' “<®n ' ^n^ “m.n-1 

(42) 

\«,n * Vl.n' * 

(43) 

'm 

(44) 


III. 5. Method of Solution. Implementation, and 
Cowputer Prooraiii 


To obtain the difference equation of eq. (4d). we 
evaluate the equation at point (l,n>i). The Inte' 
gral on the r1ght*hand side Is evaluated by using 
the trapezoidal rule. The resulting equation Is 


“l.n " ^*^l,n-l “l.n " “l.n “l,n-l 

- 4(n*3/2)juj^^,2 ‘ “i.n-* “l,n-l^^ 

« (I ^ r ] 

“l.n-1 ' l.n-1 ‘l,n-*^n>2 


III. 5. 2. Subsequent Stations 

Here. Iteration tmtst be used to obtain the 
solutions for h^^ „. „ and 

“arfl n' shape factor which signifies the 

expansion of the viscous core In the downstreaw 
direction can be determined by Iteration or by 
assuilng that the outer boundary Is a stream sur- 
face. If the shape of the outer boundary Is 
prescribed, then Is already known. Next, we 

describe the steps used to proceed In the downstream 
direction using equations (27). (29). (31) and (39). 


1. Equation (27) Is used to calculate h^^ 

Starting at n ■ 2. we march toward n > N f 1. On 
the right hand side h^^„. u^^„ 

“m.n-1 Vl.n Vl,n-1 

are unknown. The latter values sre Initially 
extrapolated from the preceding station values u, 

and u. 


m,n-r 


m.n 

(During the Iteration cycles, they are 


assmed equal to the recently calculated values of 
the preceding Iteration.) 


III. 5.1. Initial Axial Station 


Once the axial- and swirl -velocity profiles are 
specified, equation (l6), the initial pressure and 
Initial radial -velocity profiles are calculated from 
equations (10)-(13). For the Initial pressure, 
equation (11) 1$ evaluated at point (l,n-l) and we 
obtain 


'l,n-l 


P, 


.n 


404^ (n-3/2)^ 


(45) 


This equation Is used to calculate the initial 
pressure. 


For the Initial radial velocity, we eliminate 
the axial derivatives In equations (lO)-(13) end 
obtalii the following equation 


2 t 

u 3H . 3 u 1 / 1 \ / T 


(46) 


11. Equation (31) Is used next to calculate k^j^ ^ 

Following the Tridlagonal Algorithm, we obtain the 
following recurrence equations. 


*n ’1 


D - C I 
n n n-j 


n 


n n-l 


ll - 


. _ 

8n”*^7^'n-7' 


(49) 

(50) 


*'m+l,n “ ^n *^n ’Wl.n-rl 


Starting at n « 2, we march toward n » N to 
calculate I and J using eqs. (49) and (50). Since 
the value of Is known from the boundary 

conditions, we march backward from n • N to n ■ 2 
to calculate ic , using eq. (51). The unknown 
value of u^^ ^dq. (37), Is Initially extrapolated 
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frw tht prtctdins station v«1ut . (0ur1n9 the 
Ittratlon process, they ere assuntorequal to tht 
recently celculettO values of the prtceOIno Itera- 
tion. ) 

111. Equation (29) Is used next to calculate 
^•*1 n-l ***'^^"9 n • N ♦ landaarchlng toward 
n • 2. 


1v. Equation (3S) Is used next to calculate 

The following recurrence equations are used: 






-A. 


«n * 




(52) 




Lj • 0, Nj • 1 


(S3) 


Vl.n * '“n *"n Vl.n+1 


Here, we aarch upward and bactaiard In che saae way 
as step 111 using equations (52) - (54). 

v. Steps 1-1v are repeated until the calculated 
values converge. 

v1. Next, we March to the next axial station. It 
should be noted that If convergence failed at any 
or took large maaber of Iteration cycles, this 
Indicates that the breakdown point may be reached. 
The boundary-layer-llkc equations are not valid 
In this region and one has to resort to the Navler- 
Stokes equations. 


III. 5. 3. Coaputer Prograe 

A coMputer prograe is written to Impleaient the 
developed method of solution. The Nain Prograai 
reads the boundary conditions and Initial conditions 
either in the form of discrete values or In the 
fora of analytic expressions. It also reads the 
step size Oc> aspect ratio of the grid A., nuaber 

a* steps N In the ; direction, the paraawter a, 
and the ratio R. The main prograa Interacts with 
several subroutines. 

TWo subroutines are devoted to the Initial 
conditions; the first calculates the Initial 

pressure profile P, , , while the second calculates 
1 ,n-i 

the profile of initial radial velocity h, 

1 ,n 

Another five subroutines are used to calculate 

Vl. ^oH>l,n' •‘iBvl.n* ^arrl.n-l Vl.n. 

subroutines are successively called in the main 
prograM within a continuous iterative loop until 
the values of h, k. p and >j converge at two succes- 
sive cycles of Iteration. If convergence is not 
reached at any downstream station a^tar a prescri bed 
nMber of iterations, the program stops. This Is 
an indication that this station is closely upstream 
of the breakdown station. If the program converges 
at all stations, it is stopped after a prescribed 
distance in the axial direction. In the case of 
the ’eading-edge vortex core, the program stops 
after a distance equals to one and one-hald the 
wing root chord. 


A Plotting Subroutine is also included to 
plot the profiles h, k, p and u at any desired 
station. 

Figure 18 shows a flow chart of the main fea- 
tures of the computer program. 


III. 6. Nueerical Examples 

Preliminary numerical examples have been con- 
sidered to check the developed computer program. 

In this paper, the resu’ts for a tral ling-edge 
vortex core calculations are presented. Typical 
velocity profiles for the Initial viscous-core flow 
are used. This example was considered earlier by 
Hall^. The application to a leading-edge vortex 
core, with boundary conditions obtained from the 
inviscid model. Is currently considered. The OEC- 
10 Computer of the Old Cominion University C^uter 
Center is used to carry out the calculations. 

For the trai ling-edge vortex core, the follow- 
ing Initial conditions are used: 

At x^ - 0.25, we have 
u(xj.c) • 1 - 0.25 exp (-c^), 
k(x^,c) ■ 0.5[l-exp(-;^)] 

For the boundary conditions, the following 
conditions are laed: 

At ■ 6.0 and for 0.25 < x ^ l.O, we have 
u(C(j.x) • 1.0, k(c^j.x) ras, p(;^.x) • 1.0 

The shape factor of the outer boundary is 
considered constant; \(x) ■ 1.0. In the axial 
direction, 15 stations are taken (Ox ■ 0.05), 
while In the radial direction, 80 stations are 
taken (Ac ■ 0.075). The factor R is taken to be 
1.0. Figure 19 shows the axial- and swirl -velocity 
profiles at x • 0.25, 0.5, 0.75 and 1.0. One can 
see the strong Interaction betwMn the axial and 
swirl velocities. As the swirl velocity decreases 
in the downstream direction, the axial-velocity 
deficit 1-u decreases too. These results agree 
exactly with those of reference 48. The results 

of the linear theories®^ are also included in 
the figure. The computational time for this case 
took 4.0 seconds of CPU time. 

IV. Concluding Remarks 

Integrated nuMrIcal methods are developed for 
the inviscid and viscous solutions of highly swept- 
back wings at large angles of attack. The results 
of the NWV-method are remarkably successful in 
predicting Humeel's experimental data. The devel- 
oped viscous program, based on the first order 
bouncary- layer- like equations, successfully repro- 
duces Hall's nuimrlcal results. Currently, accurate 
near-field calculations of the inviscid model are 
developed to provide the boundary conditions '‘eeded 
for the viscous solutions of the leading- and 
trai ling-edge cores. 
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Figure 1. Vortex forewtlon behind a slender delta 
Min; (schematic) a) EHe and Jones 
(Ref. 8). Hummel and Redeker (Ref. 9) 
b) Huinmcl (Ref. 5). 



Viscous-Vortex Core 



b. MODEL NO. 2 



Figure 3. Typical solution of the leading-edge 
vortex sheet, AR ■ 1, 12 x 12 lattice; 
Kandll, et.al. (Ref. 36). 





Figure 2. Models of Viscous-Vortex Cores 


Figure 4. Height and spanwise oositlon of the 
leading. edge core, AS • 1, a • 15®; 
Kandll , et.al . (Ref. 36) . 
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Pigura 8. Convargad solution showing tha 1aad1ng> 
and trail Ing-adga coras In two- and 
thrat-dimnslonal vlaws; A8 « 1 , 
a - 20.5°, 12 a 12 lattica, 1/4 turn, 
NNOV-Nfthod. 



LEG TEC 

ITER • 6 




12 X 12 lattica, 1/4 turn 





12 X 12 lattica, 1/2 turn 



12 * 12 lattice, 3/4 turn 


Figura 9. Convargad solution showing the laadlng- 
anci trail Ing-edga coras In two- and 
thraa-'dlaenslonal views; AR ■ 1, a ■ 
a ■ 20. S°, 12 X 12 lattice, 1/2 turn, 
MNOV-Hati.od. 


Figure 10. Lciciti'.g- and trail ing-adjia sheets and 
flow direction in a plane perpendicular 
to tha wind direction wj(th different 
turns, AR ■ 1 , a ■ 20,5"', l • 0.08, 

IV^^l ■ 1. 
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12 X 12 littict, 3/4 turn 
Liad 1 n 9 - ind trail 1ng>adgt shafts and 
flow dlractlon In a plana parpandicular 
to tha wind dlractlon with dUfarant 
turns, AR • 1 , a • 20.5°. C ■ 0.5333. 
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12 X 12 lattice, 3/4 turn 
Flgura 12. Laadlng- and trail Ing-adga shafts and 
flow dlractlon In a plana parpandicular 
to tha wind dlractlon with diffarant 
turns AR • 1, a ■ 20.5°, t • 1,066. 
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Figurt 15. Sp»nw1$« ntt turfuce orisiurt vtrljtlon 
a - 20.5®, 12 X 12 lattice. 
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Figure 17 . Finite-Oifference grid lyitem. 
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STEADY AND UNSTEADY INCOMPRESSIBLE FREE-MAKE ANALYSIS 
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ABSTRACT 

The first part of this lecture covers detailed descriptions of the 
flows around highly sweptback wings and bodies of revolution at high 
angle of attack. Next, inviscid model approximations and mathematical 
formulation of the problem are given for steady and unsteady Incompres- 
sible flows. A general presentation of the methods of solution is 
given with emphasis on the modern computational techniques. Detailed 
descriptions of the nonlinear vortex-lattice and vortex-panel techniques 
are presented to show how the boundary conditions are enforced using 
Iteration. 

Typical numerical results are compared with the available experi- 
mental data. The lecture ends with concluding remarks regarding the 
existing computational techniques and the outstanding problems in 
this area. 
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IN RECENT YEARS, accurat« analysis and reliable prediction of vortex 
flotvs have become extremely Important to aerodvnawicists dealing with 
modem designs of fighter aircrafts, ^ missiles* and helicopters. 3 
Modern fighter aircrafts fly at high angles of attack during take-off, 
offensive and defensive maneuvering, approach and landing. In this 
range of angle of attack, vortex flows develop around tiM aircraft 
with dominant effects on Its aerodynamic characteristics and controlla- 
bility. Modern designs of missiles require high launch angles of attack 
and high maneuverability within which a very complex vortex flow 
develops. For helicopters, the Interaction of a blade with the vortex 
wake of another blade affects its operating performance, vibration and 
noise characteristics. In forward speeds, blade slap, a predominant 
source of external noise, occurs due to the rapid time rates of change 
of the blade pressure developing from Its passage through a tip vortex 
of a preceding blade. Examples of these vortex flows are given In 
Figure 1. 

For all these applications, one has to deal with strong nonlinear 
aerodynamics. Compressibility and separation of the flow are the main 
sources of the strong nonlinear effects. As the flow Mach number changes 
from low subsonic to transonic and supersonic, the flow undergoes 
several qualitative changes a.'id shc^k waves (detached or attached 
appear In the flow. As the angle of attack Increases frcxs low to 
moderate and high values, the lift- and pitching-moment coefficients 
become nonlinear functions of the angle of attack. This nonlinear 
behavior of the total loads is attrib' >d to flow separations, from 
the wing edges (side and leading edge. nd the body leeward side. In 
the fonn of vortex sheets that roll up into strong vortex cores. 

Prediction of the coupled effects of nonl Inearl ties; compressible 
and separated flow nonlinear! ties, is currently receiving considerable 
efforts from researchers working in tht computational fluid dynamics 
area, but it Is far from being complete. In the supersonic flow regime, 
a few successful attempts based on the vinite difference solution of 
Euler's equations exists. in the transonic flow regime, nothing is 
available yet and several efforts have just started. 

In the present lecture, I am going to discuss the nonlinearity due 
to flow separation only. Compressibility effects will be presented 
using the Prandtl-Glauert transformation. The full nonlinear com- 
pressibility effects will be presented in the next lecture using the 
Integral equation formulation (Green's function solutions}. In this 
regard, the lecturer has presented'O a steady nonlin^iir vortex lattice 
technique which accounts for the full nonlinear compressibility terras 
In the nonlinear potential equation. Flow compressibility in this 
technique is modeled by volume-source distribution within a finite 
volume around the body. The technique has not been numerically tested 
yet but work is underway using a similar technique which includes 
shock capturing using the artificial compressibility meth<Ml of reference 
11 . 

The first part of this lecture covers a description of the flow 
based on the available experimental measurements. This is followed by 
the inviscid model approximations and the mathematical formulation of 
the steady and unsteady problems. A general presentation of the methods 
of solutions is given with emphasis on the modern computational techni 
ques. Detailed descriptions of the nonlinear vortex-lattice and vortex- 
panel techniques are presented to show how the boundary conditions are 
enforced using iteration. Typical numerical results are compared with 
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the available experimental data. The lecture ends with concluding 
remarks regarding the existing computational techniques and the out- 
standing problems in this area. 

Major portions of this lecture are based on earlier published 
papers* reports and theses 2-22 by the lecturer, his collaborators 
and his students. 

DESCRIPTION OF THE FLOW FIELD 

In the low-subsonic to transonic speed range, numerous experimental 
data are available for vortex flows about highly swept back wings23-38 
and bodies of revolution^9-51 at an angle of attack. 

HIGHLY SWEPT BACK WINGS- The fluid flow past such wings is charac- 
teri by flow separations from side and leading edges due to strong 

crosi lows The flows from the pressure and suction sides of the wing 
leave at tht side and leading edges to form free-shear layers. The 
separated free-shear layers roll up spirally and form two vortex cotes 
which are continuously fed by vorticity shed from the attached boundary 
layers on the wing surfaces through the free-shear layers. This type 
of flow separation is known as the "primary separation." It has a 
dominant effect on the aerodynamic characteristics due to the large 
strength of its vortex core and its prximity to the upper surface 
of the wing. It generates a large suction pressure peak on the upper 
surface under the primary vortex core. 

The primary separation lines nf wings with sharp leading and side 
edges are fixed along these edges and hence they are known a priori. 

In the case of rounded edges, the separation lines are inboard of the 
edges and are not known a priori. Hence a viscous model or perhaps an 
empirical formula is needed to locate these lines. 

The adverse pressure gradient outboard of the suction peaks affects 
the boundary layer flow on the upper surface and "secondary separation" 
from the wing surface occurs. The secondary separated flow forms either 
an additional free shear layer or a bubble depending on the angle of 
attack. In the range of moderate to large angles of attack, the secon- 
dary free-shear layer rolls up spirally in an opposite sense to that of 
the primary shear layer and forms a secondary vortex core with a 
strength much smaller than that of the primary core and of opposite 
sense. The effect of secondary separation depends on the type of 
boundary- layer flow on the upper surface. For a laminar boundary layer, 
the secondary core produces another lower pressure peak between the 
secondary and primary lines of separation. For a turbulent boundary 
layer, such a pressure peak is hardly noticed, and the primary pressure 
peak is slightly higher than that of the laminar boundary-! ayer flow.^ 

A third type of flow involving a "tertiary" separation may occur 
between the lines of secondary and primary separation due to the adverse 
pressure gradient created by the secondary vortex core. Figure 2 gives 
a schematic of these flow separations. 

The free-shear layer emanating from the trailing edge is of oppo- 
site strength to that of the primary free-shear layer. Within a short 
distance behind the trailing edge (a distance of a 1/4 root-chord 
length for a delta wing of aspect ratio of one at 20.5^ angle of attack), 
the trailing free-shear layer rolls up spirally in an opposite sense to 
that of the primary free-shear layer and forms two counter rotating 
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vortex cores. The trail ir.g>edge vortex core has the same sense of 
rotation and almost the same spanwise location as those of the secondary 
vortex core, although each originates from a different phenomenon. 

The former is due to an edge separation while the latter is due to a 
surface separation. Hummel4 concluded from his latest experimental 
measurements for turbulent boundary layers at the trailing edge that 
the secondary- vortex core decays rapidly behind the trailing edge 
and its remains roll up into the trailing vortex. 

Predicting the correct formation of the trailing-edge vortex core 
is essential for predicting the correct aerodynamic loads near the 
wing trailing edge due to its proximity to this dege. Moreover, a 
correct prediction is of paramount Importance to problems involving 
high-angle-of-attack aerodynamics for the canard-wing configurations. 

The size of the primary-vortex core and viscosity of the flow with- 
in a narrow region around the core centerline play important roles in 
the aerodynamic characteristics of low-aspect ratio wings at large 
angles of attack. For these types of wings, it is well known that the 
lift coefficient increases nonlinearly with the increase in the angle 
of attacM. This is attributed to the increase of the strength of the 
primary vortex core. However, the increase in the lift coefficient 
is limited by the occurrence of vortex-core breakdown over the wing and 
hence there exists a maximum angle of attack corresponding to the maxi- 
mum lift. Figure 3. 

The phenomenon of vortex-core breakdown is defined as an abrupt 
increase in the core diameter. It occurs due to the adverse swirl - 
induced pressure gradients along the core axis. Using the radial com- 
ponent of the momentum equation of the boundary-layer like equations 
(for quasi -cylinderical vortex cores), one can show that an increase 
in the swirl velocity of the flow, due to an increase in the core 
strength, int.'eases the axial pressure gradient along the core axis. 
Therefore, the flow near the core axis loses its axial momentum and 
swelling of the vortex core develops. 

For moderate angles of attack, vortex-core breakdown develops far 
downstream from the wing trailing edge. As the angle of attack 
increases, vortex-core strength increases and the breakdown point 
travels upstream toward the trailing edge. As the breakdown point 
approaches the trailing edge, the slope of the lift versus angle of 
attack curve decreases until the maximum lift' is reached. Thereafter, 
the lift coefficient drops sharply, the breakdown point crosses the 
trailing edge and the wing stalls. 

SLENDER BODIES-Three distinct patterns of flow separation develop 
on the leeward side of a slender body as it is pitched through three 
ranges of aagles of attack. 

At low to moderate angles of attack (6°-25®), two symmetric vor- 
tices appear on the leeward side which consist of two counter rotating 
vortex cores with increasing strength and size in the downstream direc- 
tion. Vorticity from the boundary-! ayer flow on the body surface is 
shed from the separation lines and feeds the vortex cores through roll- 
up free-shear layers. This type of flow separation is also known as the 
primary separation. The separated flow has a dominant effect on the 
pressure distribution on the leeward side. Similar to the flow around • 
delta wings, secondary-flow separation develops also inboard of the 
primary separation line. 
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Above the 25® angle of attack, the two continuous feeding free- 
shear layers break>up asynmetrical 1y at points along the two separation 
lines and multiple pairs of asymmetric vortex cores form over the lee- 
ward side of the body. Each vortex core originates from a separate 
point on the separation line. Ericsson and Reding52 conclude that both 
the "vortex crowding" near the body nose and the degree of freedom of 
the separation point may contribute to the formation of asymmetric 
vortices. The resulting asymmetric pattern of vortex cores produces 
large side forces which can be larger than the normal force. 

At very large angles of attack, above 60® angle of attack, the 
crossflow becomes dominant and unsteady vortex shedding starts. The 
separated free-shear layers are shed in the form of a Karman vortex 
street. 

APPROXIMATIONS FOR INVISCID ANALYSIS 

The experimental results given above show that we are dealing with 
highly complex flow fields. For highly swept back wings, four different 
vortex cores {primary, secondary, tertiary, and trai ling-edge) have been 
found and apparently their origins are well understood. For bodies of 
revolutions, two different vortex cores (primary and secondary) have 
been found and three different modes of flow separation (steady sym- 
metric, steady asymmetric and unsteady asyninetric) have been identified. 
For inviscid analysis, certain approximations are adopted so that the 
resulting potential flow model represents the main features of the real 
flow to a good degree of accuracy. 

For wings, the effect of the tertiary vortex core is very small and 
can be neglected. The effect of the secondary vortex core is small in 
symmetric flows particulary when the boundary layer on the suction side 
is turbulent. Hence, its effect on surface pressures and integrated 
forces can be neglected without appreciable error under the conditions 
stated. Thus, we are left with the primary flow separation and the 
trailing-edge flow separation. For an isolated wing, the effect of 
deformation of the primary and trailing-edge vertex cores and their 
feeding free-shear layers is neglected after a distance of about 1/2 
to 1 root chord since it is a far field effect. For canard-wing, 
wing-tail and strake-wing configurations, this effect must be taken 
into account. 

In most inviscid analyses, the attached boundary layers are 
replaced with bound-vortex sheets while the free-shear layers are 
replaced with free-vortex sheets. Accordingly, the inviscid region 
representing the vortex core should be represented by a conical sprial 
vortex sheet which ends with a concentrated vortex line along the core 
centerline. Although this theoretical representation of the vortex- 
core is ideal for the inviscid model, the usage of many turns within 
the core region creates numerical problems in a three-dimensional model 
which does not assume slender-body approximations. The numerical 
problems arise due to the proximity of the vortex layers of the spiral 
and due to the large number of small -size panels needed to accurately 
model the turns. However, it was found, according to the numerical 
results and their excellent comparison with the experimental data, 
that only one turn of the spiral need to be accurately modeled while 
the remaining turns of the spiral are replaced jy a cut ending with a 
concentrated vortex line along its edge (Figure 4.) 
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Since one does not know a priori where the tralllng-edge vortex core 
originates, one has to obtain it as a part of the solution. This has 
been successfully accomplished, for the first time, by the lecturer '3,15. 

In all inviscid wing models, the wing edges are assumed to be 
sharp so that the separation lines are known a priori and hence viscous 
modeling is not needed to determine the lines of separations. However, 
for wings with round edges, viscous modeling is needed to determine 
these lines. Moreover, vortex-breakdown points are assumed to be far 
downstream so that the variations in the size of the primary-vortex 
core can be neglected. In fact, this assumption limits the large 
angles of attack at which -inviscid modeling is applicable. With the 
vorticity confined to the inner boundary of the flow region (bound- 
and free-vortex sheets), the resulting model is a potential-flow model 
governed by Laplace's equation (in i incompressible flows) and certain 
boundary conditions. 

For bodies, the effect of secondary separation is also small and 
can be neglected. Unlike the primary separation of wings, the primary 
separation lines are not known a priori and hence viscous solutions 
are needed to locate these lines nd i'»termine the strength of shed 
vorticity (this is a separation from a smooth surface and hence Kutta 
condition is not applicable). In this regard, the reader is referred 
to the work of Mendenhall, Spangler and Perkins. 53 Since we are 
dealing with inviscid analyses only, this problem will not be discussed 
anymore in the present lecture. The rest of the lecture addresses the 
steady and unsteady problems of wings with vortex flows. 

MATHEMATICAL FORMULATION OF THE GENERAL PROBLEM 

The problem is formulated relative to a wing-fixed frame of refer- 
ence xyz. The x-axis is the wing centerline and the xy-plane is the 
wing plane of symmetry. The wingis rotating at the angular velocity 
S and the freestream velocity is 0„. The general orientation of the 
wing is described using the Eulerian angles a, e, and y» Figure 5, which 
refer to the angles of attack, yaw and roll; respectively. In terms 
of these angles and their time rate of change, the dimensionless 
freestream velocity and the angular velocity are expressed by 

e^ * cosa COS6 T + (-sincx cosy ♦ cosa sing siny)3 + (sina siny 

♦ COSO sina cosy)ic (1) 

n * (-0 sing + y)i + (a cosa siny + g COSy)j + (a cosg cosy 

- 8 siny)k = + n^k (2) 

where i, j and x are the base unit vectors of the xyz-frame of 
reference. 

The unsteady irrotational ideal flow in the region outside of the 
wing and its separated vortex sheets is governed by the Laplace's 
equation 
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■ 0 (3) 

where <^{r,t) is the perturbation velocity potential. The no-penetra- 
tion condition on the wing s(r,t), relative to the wing fixed frame 
of reference, is given by 

+ (e + V(ji - nxr)*vs » 0 on s(r,t) ■ 0 (4-a) 

For a rigid wing, ||- » 0 and equation (4-a) reduces to 

(i + 7(|i - Sxr)*n =0 on s(r) * 0 (4-b) 

* s 

On the separated free-vortex sheet w(r,t), the no-penetration condition 
is given by 

Iy + " Sxr)*7w » 0 on w(r,t) * 0 (5) 

The no-pressure jump condition on w(r,t) is obtained from the unsteady 
Bernoulli's equation 

Cp(r.t) = -7^.[7(j» + 2(e^ - axr)] - 2 || (a) 

where C (r,t) is the pressure coefficient at any point r and at any 
time t.^ Forming the pressure jump from equation (a) and equating the 
result to zero, we obtain 

iC ■ C , - Cp 2 - + 2(i, - 5xr)] 

- 2 |i*l ■ *2> • 0 (b) 

where the subscripts 1 and 2 refer to the upper and lower surface of 
the wing; respectively. Rearranging equation (b) and setting 

<|i^ - <j >2 » A<j> (c) 

one obtains 

where v^ * velocity of a wake element relative to the xyz-frame of 
reference 

- i ( 7 ^^ + 7.^2) + K ■ 


Equation (d) represents the theorems of Klelvin and Helmhotlz 
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5T * §t * ”a ^ ° O'’ w(r,t) - 0 (6) 

of conservation of the circulation and the outflow of vorticity; 
respectively. In equation (6), n^ is a unit normal to the surface A 
bounded by a closed curve around which the circulation r is calculated. 
Equation (6) simply states that the rate of change of circulation around 
a closed curve or the rate of change of outflow of vorticity through the 
surface bounded by this closed curve is zero (following the same fluid 
particles. ; 

For uniqueness of the solution, one has to impose the Kutta condi« 
tion along the edges of separation. Here. Kutta condition is repre- 
sented by 


A Cp I - 0 (7) 

te.le.se 

Finally, the infinity condition requires that 

-► 0 away from s and w (8) 

Equations(l) - (8) are the required equations for the general unsteady 
flow problem. 

For steady and unsteady symmetric flows, equations (l)-(8) are 
simplified as follows: 

STEADY-SYMMETRIC FLOW- For this problem, we set 0*0, b=y* 0 and 


drop the time dependent terms. The resulting equations are 

■ cos o ? - sin a 3 (9) 

7^* - 0 (10) 

(®, + 7(ji) • n^ » 0 on s(r) » 0 (11) 

(i ,+ V(^) • n^ » 0 on w(f) « 0 (12) 

ACp ■ -2(n^xiL) • (e^ + 7<|»)«0 on w(r) » 0 (13) 

ACp I *0 (14) 

te.le.se 

7<^ 0 away from s and w (15) 

UNSTEADY-S^^EITIIC FLOW-For this problem, we set 6 *y*B*y* 0 and 
the resulting equations are 

• cos o ? - sin a j (16) 

5 - 0^ ic » a k (17) 
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• 0 



08) 

(e. 

+ . 

«,xr) • ", 

■ 0 on s(r) ■ 0 


(19) 

3t 

(e + Vi(i - n,xr). • 
• 2 

7w ■ 0 on w(r,t) ■ 

0 

(20) 


DT . 
fft 

0 

on w(r,t) ■ 

0 

(21) 


ACp 

1 * ^ 



(22) 


TE. 

LE,SE 





7^ -► 

0 

away from s and 

w 

(23) 


EXISTING METHODS OF SOLUTION 

Before the Nonlinear Vortex-Lattice [Nonlinear Discrete-Vortex (NDV)] 
and Nonlinear Vortex-Panel [Nonlinear Hybrid-Vortex (NHV)] methods are 
discussed in detail, a brief discussion of the existing mathematical 
models is presented. The literature contains several steady and unsteady 
inviscid-flow models with various degrees of limitations and drawbacks. 
These models can be divided into four main groups. Each group is 
presented and evaluated in the next subsections. 

SLENDER-BODY MODELS-This group of models uses slender body and 
conical flow assumptions. 54-61 Modeling of the primary-vortex core and 
its feeding sheet was first introduced by Brown and Michael. 55 However, 
the feeding sheet in their model was taken as a planar surface and hence 
it did not represent the real rolled-up vortex sheet. Mangier and 
Smith56 introduced the first realistic model of the primary-vortex core 
and its feeding sheet. However, this model does not account for the 
feeding-sheet deformation in the chordwise direction due to the slender- 
body assumption. These models satisfactorily predict the pressure 
distribution on the front portion of the wing surface. In the rear 
portion, the models fail to predict satisfactory pressure distributions 
because kutta condition cannot be satisfied at the trailing edge. Such 
models were limited to slender delta pi anforms. 

OLD NDV-MODELS-In these models ,'0*17-22,62-66 the bound-vortex 
sheet and the free-vortex sheets are approximated by a set of concen- 
trated vortex lines. The bound-vortex sheet is replaced by a bound- 
vortex lattice, while the free-vortex sheet is replaced by segmented 
free-vortex lines (in the case of steady flow) or by a growing free- 
vortex lattice (in the Las? of unsteady flow). The boundary conditions 
are satisfied at certain control points on the bound- and free-vortex 
system using an iterative technique. Excellent agreement was foundl0*22 
between calculated and experiri^ental total aerodynamic characteristics, 
and the agreement between calculated and experimental total aerodynamic 
characteristics was satisfactory for wings with only side-edge separa- 
tion. For wings with leading-edge separation, however, the agreement 
was less than satisfactory for some cases. Although the discrete- 
vortex model has for many years worked very well for attached-flow 
problems, 67 when vortex-type separation from leading edges and/or tips 
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occurs, the free-vortex system lies close to the lifting-surface and 
the results are found to be sensitive to the model of the primary 
separated flow. The most obvious drawback of the old model Is the 
lack of a realistic model of the primary vortex core and Its feeding 
vortex sheet. 

Although the old NOV-method was almost abandoned In predicting the 
flow details and the distributed aerodynamic characteristics, the 
modified NDV-method^3,15 pinpointed and cured most of the problems 
encountered with the old method. It enjoyed a remarkable success In 
predicting, for the first time, the latest experimental results of 
Hummel co 7c 

DOUBLET-PANEL MODELS-Thls group of models employs doublet-panels”**" 

In this TOthod, the wing and Its free-vortex sheets are divided into 
networks of quadrilateral panels. Each panel of the networks repre- 
senting the wing has a biquadratic local doublet distributions and a 
bilinear local source. distribution. The panels of networks representing 
the free-vortex sheets have biquadratic local doublet distributions. 

Source and doublet splines are used to express the distributions of 
singularities on the networks in terms of discrete values of singularity 
strength at certain standard points on each network. The boundary 
conditions and continuity of singularity strengths across abutting 
networks are enforced at certain standard points on network. The 
results of this method’3,/4 generally good when the solution con- 
verges. Apparently, the difficulty in obtaining convergence is due to 
the failure in satisfying the continuity of the derivatives of the 
doublet strength across aL-ittlng networks. This is equivalent to the 
existence of concentrated vortex lines between abutting networks. 

The doublet panel method was extended^S to include the effect of 
entrainment of the primary vortex cores through an empirical approach. 

The results indicated that the entrainment increased the normal-force 
coefficient substantially over the experimental values. 

This group of models do not account for the chordwise deformation 
of the primary-vortex sheet. Consequently, the primary vortex sheet 
cannot be fed three-dimensionally with the shed vorticity. Moreover, 
the trailing-edge free-vortex sheet was represented by a "fixed design 
wake." 

NHV-MODELS-Thls group of models employs a nonlinear hybrid vortex 
method12,14,16. this method, continuous-vorticity and vortex-line 
representations of the wing and its separated free-shear layers are 
used. Continuous vorticity is used in the near-field calculations 
while discrete vortex-lines are used in the far-field calculations. 

The winc/ and its free-shear layers are divided into planar quadri- 
lateral panels having first-order vorticity distribution. The aero- 
dynamic boundary conditions and continuity of the vorticity distributions 
are imposed at certain nodal points on the panels. To satisfy these 
conditions, a? iterative technique is followed which alternatively yields 
the local vorticity distribution on the bound-vortex panels and the 
shape of the free-vortex panels. This method has been used to calculate 
the steady distributed and total loads on planar-low-aspect-raio rectan- 
gular wings. The results have shown that the spanwise variations of the 
load coefficients are in good agreement with the experimental data. 
Comparisons of the results with those of the NDV -method have shown that 
the hybrid method requires less number of vortex panels for the same 
accuracy. 
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NONLINEAR OISCRETE'VORTEX (NDV] METHOD 

In this method the basic unknowns are the circulation distribution 
and the shape of the free-vortex sheets. The following sub-sections 
give the details of the method. 

CONSTRUCTION OF THE DISCRETE- VORTEX METHOD-FIgure 6 shows how the 
discrete vortex model Is constructed for a delta wing. Although the 
example discussed here Is for a thin, flat, delta wing, the method Is 
general and Is not restricted by the geometrical parameters of the wing; 
e.g. camber, aspect or thickness ratios or wing planform. 

The first fitep Is to div1;1e the wing into rectangular and cropped- 
delta winglets as shown by the dashed lines in Figure 6. a. A rectan- 
gular winglet is aerodynamical ly represented by a spanwise bound- 
vortex segnent of constant circulation ri . This segment is placed at 
the quarter-chord length of the winglet (the chord length of the 
rectangular winglet is the characteristic length of the problem). 

In addition a control point is placed at the three quarter-chord length. 
The choice of these positions is suggested by thin airfoil theory. 

It can be shown that the bound-vortex sheet representing the two- 
dimensional flow around a f:at plate at an angle of attack can be 
replaced by a point vortex of the same strength as that of the continuous 
vortex sheet under the following conditions: a) the point vortex is 
placed at the quarter-chord length and b) the flow tangency condition 
Is enforced at only one point at the three-quarter-chord length. 

On the other hand, a cropped-delta winglet is aerodynamical ly 
represented by a bound-vortex segment of constant circulation. This 
vortex segment is directed along the perpendicular from the midpoint 
of the winglet root chord to its leading edge. With this choice it can 
be seen that the vorticity of this vortex segment does not have a 
component along the leading edge and hence the Kutta condition is 
approximately satisfied along this edge. 

Chordwise bound-vortex segments arise due to the differences in 
the strengths of the neighboring spanwise, bound-vortex segments. In 
this way, a bound-vortex lattice which replaces the continuous, bound- 
vortex sheet is constructed. The model is completed by adding free- 
vortex lines, representing the continuous free-vortex sheets at the 
ends of the bound-vortex lattice along the edges of separation > the 
leading and trailing edges. Each line is divided into a series of 
small, straight segments (near-wake region) and one semi-infinite 
vortex line (far-wake region). The upstream end of each segment 
represents a control point of the wake surface where the kinematic and 
dynamic boundary conditions are satisfied. The initial positions and 
shapes of these lines are prescribed. The resulting model is shown in 
Figure 6b. This model has an unknown circulation distribution and a 
wake that can be deformed to satisfy the boundary conditions. 

The model described above is used to solve the steady-flow problem 
by satisfying the corresponding boundary conditions. On the other hand, 
If the problem under consideration Is for an unsteady flow which starts 
from a steady flow situation, then the solution of the steady-flow 
problem serves as an initial condition to the unsteady problem. Further- 
more, if the problem under consideration is for an'unsteady flow which 
arises from an impulsive motion of the wing, then the irtlal condition 
corresponds also to the solution of the model given above, but with the 
wakes renx>ved from the model. 



12 


CALCULATION OF THE VELOCITY FIELD-To satisfy the boundary conditions 
on the wing and Its wake and to calculate the surface pressure distri- 
bution, one needs an accurate method to calculate the velocity at any 
field point rj at any time step tk. If the field point is off the wing 
and Its wake, then the velocity Is given by 

V(Fj,t,^) “ 7(^(rj,tj^) + - n(t^) X (24) 

where 

n(tk) 

^ [r^(t|j)/4irh^(rj,t^)][cos0i^(rj,tj^) 

- cos02^(rj,t^)]i^(rj,t|^) (25) 

Is the Induced velocity from a11 the vortex segments of the model. The 
parameters on the right-hand side of equation (25) are those of Blot- 
Savart's law^ff. The number of vortex segments n(tk) is a function of 
the time step tk due to the growing vortex lattice In the wake In the 
unsteady-flow problem. To account for the near-vortex velocity, equa- 
tion (25) Is modified by an artificial viscosity which Is obtained 
from reference 77, 

When the field point is on the wing surface or on the wake surface, 
one has to account for the self-induced tangential velocity due to the 
local strength of the vortex sheet. Figure 7 shows the parameters 
Involved in calculating the components of the induced tangential velo- 
city In the x and z direction at a point p for a rectangular vortex 
element In the xz-plane. With linear interpolation. It Is easy to ^how 
that these components are given by 

Vtx(y • 0-) * - 0/2 tj)[ri(x + - Xi) + F3 (xi - x)]ei (26. a) 

Vt^Cy • 2’^ ‘ ^ 0/2 4)[r4(z + iz - 2i) + r2(zi - z)]e2 (26. b) 

where the arguments y » and y » O" correspond to the upper surface 
and the lower surface of the wing, respectively. Equations (26) must 
be added to equation (24) if oneis to calculate the pressure distri- 
bution on the upper and lower surfaces by using Bernoulli's equation. 
Extension of equations (26) to a general, quadrilateral vortex element 
is straightforward. 

BOUNDARY CONDITIONS FOR THE STEADY PROBLEM- Jhe boundary conditions 
on the wing surface s(r) and the wake surface w(r) are satisfied by an 
Iterative process. To initiate the iterative process, one needs to 
prescribe an initial geometry of the wake surface. It has been found 
from several numerical tests that the number of iterative cycles 
required to achieve the solution can be reduced by an appropriate choice 
of *'he initial geometry. This initial geometry depends on the problem 
under consideration and thus it varies from one problem to the other. 

For Instance, the number of iterative cycles for the steady, 
symmetric-flow problem is reduced by about 20% when the free-vortex 
lines emanating from the leading edge are prescribed to be straight 
lines pitched at one half the wing angle of attack. In addition, 
those lines emanating from the trailing edge are assumed straight lines 
pitched at one third the wing angle of attack. Here, the comparison is 
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made with respect to the number of Iterative cycles required for the 
same problem when all the free-vortex lines are prescribed to be 
straight lines pitched at an angle equal to the wing angle of attack. 

In the case of a steadily, rolling wing at zero angle of attack, 
an approprlate^lnitlal guess is found to be related to an angle e(r} • 

± 1/2 tan"^ |nxr|/U«. Here, we specify the free-vortex lines emanating 
from the edges c.f the advancing and receding sides t be straight lines 
pitched at the angles +e and - 0 ; respectively. 

Next, the flow-tangency condition and the spatial conservation of 
circulation are satisfied at the control points and node points, 
respectively, of the bound-vortex lattice. Thus, we obtain a set of 
linear algebric equations which yellds the circulation distribution r.j. 

With the circulation distribution fixed, the kinematic and dynamic 
boundary conditions at the control points of the free-vortex lines are 
satisfied. For steady flows, these two conditions are combined into a 
simple condition in which we require that each vortex segment In the 
wake be aligned with the local velocity at its upstream end (a control 
point on the wake surface). This means that each vortex segment Is a 
segment of a streamline (kinematic condition). Moreover, It means 
that the force on each vortex segment Is zero according to Kutta- 
Jawkowskl theorem in the small (dynamic condition). This process Is 
carried out by calculating the downstream end of each vortex segment 
according to 

where rj and rj+i are the position vectors of the_upstream and downstream 
ends, respectively, M Is the segment length and Vj is the velocity at 
Its point [equation (z4) for steady flowsL 

The Iteration sci.eme moves back and forth from the control points 
of the bound-vortex lattice to the control points of the free-vortex 
lines until convergence is achieved. We consider the iteration scheme 
converged when the variation In the circulation distribution or the 
displacement of the downstream ends of the free-vortex sepent between 
two seccessive Iteration cycles does not exceed a certain prescribed 
tolerance. Once convergence Is achieved, we calculate the pressure 
distribution and the total load coefficients. 

BOUNDARY CONDITIONS FOR THE UNSTEADY PROBLEM-Here , we consider the 
problem of unsteady flow which starts from a steady flow situation. The 
continuous motion of the wing is discretized into a series of Impulsive 
changes occuring at discrete time steps. At each time step t|(, a set 
of starting vortices develops along the edges of separation and are shed 
with the local velocities to restore the smoothness of flow at the edges 
(Kutta condition). In the same time, the starting vortices shed in the 
wake at earlier time steps are convected downstream with the local 
velocities without changing their strengths. This process satisfies 
the kinematic condition on the wake (a wake element moves along the 
direction of the local velocity) and it also satisfies the dynamic 
condition on the wake (a wake element satisfies Kelvin-Helmhol tz 
theorem). 

The position of any shed vortex rj at any time step t)( is determined 


by 
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Mhtrc tj^.i is the preceding time step and Vi Is given by equation (24). 
The strength of any newly shed vortex Is related to the change In the 
bound circulation. Hence, with the positions of the shed vortices 
known from equation (28) and with the strength of the newly shed 
vortices given in terms of the change of the bound circulation, the 
flow tangency condition at the control points of the wing yields the 
unknown circulation distribution. Tg account for the error in equation 
(2S) (because of using the velocity Vj(ri,t|(_i) at the preceding time 
step t)(>] rather than the current time step t’<), an iteration procedure 
similar to that of the steady-flow problem is performed. 

In both the steady and unsteady flows, the only difference between 
the symmetric and asynsnetric problems is the longer computational time 
required for the latter problem as compared to that of the former pro- 
blem. In the former problem, we need only to use half the wing to 
Obtain the solution because of the flow symmetry. In the latter problem, 
the whole wing must be used to obtain the solution. 

CALCUUTION OF THE PRESSURE COEFFICIENT- The distribution of the 
pressure coefficient on the upper and lower surfaces of the wing is 
calculated by using Bernoulli's equation in terms of a wing-fixed frame 
of reference. 




- [V(fj,t^)]^ + 2V(?j,t^) • [n(t^) X fj - e^J 

34(r.,t. ) 

2 J ^ 


where 


5(r^.V • ^ 


(29) 

(30) 


rj is the position vector of the control point, the positive and nega- 
tive superscripts refer to the upper_surface and lower surface of the 
wing, respectively, and v<(i, vtx «nd v^z are given by equations (25) , 

(26. a) and (26. b), respectively. The pressure is calculated at the 
control points of the bound-vortex lattice because these are the points 
where the flow tangency condition is enforced. 

In the steady-flow problem, the last term of the right-hand side 
of equation (29) is zero and all the other terms nre time Independent. 

The net presrure coefficient is given by 




- 4[v^(r],t^) t v^^(r],t^)] • [n(t^) x 


e- 


- v^(rj,tj^)] - 2[r(Fj,tj^) - (31) 

The total-load coefficients are obtained by integrating the net pressure 
coefficient on the wing. 

NUMERICAL EXAMPLES; NDV-METH0D-F1 gures 3-15 Show typical solutions 
of free-vortex sheets emanating from side, leading and trailing edges 
of isolated wings, a canard-wing configuration and a wing-body combina- 
tion for steady symmetric and asymmetric flows and for unsteady syimie- 
tric flows. 


In figure 8, the arrangement of the boumJ- vortex lattice used for 
rectangular wings Is also shown. It is noticed that the tip vortex 
covers about 13% of the wing semi -span and hence its effect on the 
spanwise pressure distribution is confined in this region. 

In figure 3, Rehbach's solutlon^^ for a delta wing of aspect 
ratio of 1 1s shown. Starting with the solution of a rectangular 
wing, he progressive! y deform the wing leading edge to obtain the 
solution of the delta wing; a time consuming process. The present NDV- 
fwsthod solves the problem directly without going through this deforma- 
tion process. Although the arrangement of the bound-vortex lattice 
of the present method differs from Rebach‘s (compare Figure 6.b and 
9.4), the final shapes of the leading-edge vortex lines are the same. 

Figure 10 shows the solution for a delta wing of aspect ratio of 
0.7. Figure 11 shows the solytion for the same wing which Is now 
steadily rolled at zero angle of attack. Notice the antisymmetric roll 
up of the leading-edge vortex lines. 

Figure 12 shows the solution for a yawed delta wing and the 
relative sizes of the leading-edge vortex sheets emanating frtjm the 
windward and leeward sides. 

In Figure 13, it is seen that the vortex trail of the canard is 
deflected Inboard of the main wing toward its plane of synmetry. It 
also moves downwards toward the surface of the main wing. This pro- 
duces a downwash on the central part of the main wing and an upwash 
on its tip region. 

In Figure 14, the solution of the leading-edge vortex for a wing- 
body combination without leeward-side separation from the body is 
1)1 ven. 

Figure 15 shows the development of the leading- and trai ling-edge 
vortex system for a delta wing undergoing a sinusoidal pitching motion. 

Figures 16-21 show samples of the computed total and distributed 
loads for the cases considered above. 

It is seen from figure 16 that the section-lift peak near the 
wing tip moves inboards as the angle of attack increases, in the linear 
case, indicated by the dotted lines, the tip vortex is not taken into 
account and one can see the substantial effect when the tip vortex Is 
included {indicated fay the solid lines). 

Figures 17 and IS show the spanwise variation of the pressure 
coefficient. The suction peak of the pressure coefficient exists 
almost at the same spanwise location of the center of the leading edge 
vortex. 

Figure IS shows the pitching-moment coefficient versus the lift 
coefficient and the lift coefficient versus the angle of attack for 
a wing-body combination. Although the leeward-side separation from 
the body is net taken into account, the computed results are in 
good agreement with the experimental data. 

Figure 20 shows the well-known hysteritic behavior encountered In 
unsteady flows as the angle of attack increases to a maximum, then 
decreases to its initial value. 

Figure 21 shows the variation of the normal-force and pitching- 
moment coefficients with time for the delta wing considered in Figure 15. 
Notice the phase lag between the loads and the wing motion. 
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THE MODIFIED NONLINEAR -DISCRETE VORTEX (MNOV) -METHOD; LEADING- AND 
TRAILING-EDGE VORTEX CORES (LEC & TEC) 

Figure 22 shows a recent converged solution of the system of free- 
vortex lines with a long deformed wake. It can be seen that the trace 
of the tralllng-edge vortex system In cross planes indicates that the 
sheet tends to deform upwards showing. a tendency to form a tralling- 
edge vortex core. However, the cross-flow planes taken further down- 
stream show that the free-vortex lines leapfrog. This does not repre- 
sent the real flow. Figure 23 shows four of these cross-flow planes 
taken perpendicular to the wind direction. 

It 1$ Cleary seen from the example given above, that the existing 
model of the NDV -method does r.>t realistically model the leading- and 
tralllng-edge vortex cores. Therefore, the model and the mxnerical 
technique must be modified In order to obtain realistic vortex-core 
•modeling. 

In the MNDV -method a realistic i. odeling of the primary-vortex core 
region Is Introduced. The old NDV-method 1$ only used In the first 
Iterative cycle to initiate the roll-up process ano to calculate the 
centroid of the leading-edge vortex system. Next, the leading- and 
tralllng-edge vortex segments are replaced by smaller segments. Then, 
the Iterative cycles proceed. 

In a typical Iterative cycle, each free-vortex line of the leading- 
edge vortex system Is allowed, three-dimensional ly, to rotate around the 
most recently calculated centroid a prescribed portion of a turn 
(1/4, 1/2, 3/4, or 1 turn). This Is done by continuously checking 
the coordinates of the free-vortex segments (as they are adjusted) 
with the location of the calculated centroid. Once this Is achieved, 
the remaining free-vortex se^ents of each line are dumped Into the 
calculated centroldal line of the vortex system. The Iterative 
technique Is followed until the circulation distribution converges. 

Figure 24 shews typical solutions at different Iterative cycles. 

It can be seen , after one Iterative cycle, that the system of free- 
vortex lines shows good ro1 1-up. The converged solution, indicated by ITER-6, 
shows the leading-edge core and Its feeding free-vortex lines. It can 
also be seen that the free-vortex lines continue to feed the LEC 
beyond the trailing edge. The tralllng-edge core Is also Indicated 
on the figure. This will be clearly seen In the cross-flow planes 
discussed In the next sub-section. 

NUMERICAL EXAMPLES-F1 gures 25-27 show comparisons between the 
experimental and numerical results for leading- and tralllng-edge 
vortex sheets and flow direction in cross-flow planes perpendicular to 
the wind direction. The numerical results are drawn at the same scale 
as that of the experimental results. The predicted sizes and locations 
of the leading- and trail ing-edge vortex sheets are In excellent agree- 
ment with the measured ones. The results of the 3/4 turn show 
the correct roll-up of the trail ing-edge vortex sheet and the correct 
locations of the LEC and TEC. This Is exactly what we expected when 
we Increased the roll-up from 1/2 to .^'4 turn. With the 3/4 turn of 
the leading-edge vertex system, the roll-up tightens and larger velo- 
cities are Induced at the tralllng-edge vortex sheet which cause It to 
deform upwards and leftwards. 


% 
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Figure 28 shows comparisons between the predicted and measured 
static-pressure contours at c = 1.066. The predicted sizes, locations, 
and levels of the pressure contours are in good agreement with those of 
the measured data. 

Figure 29 shows comparisons of the predicted and measured span- 
wise net surface-pressure variation at difference chordwise stations. 

It is seen that the results of the 3/4 turn are more accurate than that 
of the 1/2 turn. With 3/4 turn, the roll-up of the leadingedge vortex 
system tightens and produces better predicted distribution than that of 
the 1/2 turn. 

The developed computer program of the MNDV-method is computationally 
efficient. On a CYBER-175 machine, a typical solution using a 12x12 
lattice and including the cross-flow planes calculation takes 10 minutes 
of CPU time. 

THE NONLINEAR HYBRID-VORTEX (NHV-METHOD) 

In this method, vortex panels with first-order vorticity distribu- 
tion Is used in the near-field calculations. In the far field calcula- 
tions, the distributed vorticity over each far-field panel is lumped 
into equivalent concentrated vortex lines. In this way, accuracy is 
satisfied in the near field while computational efficiency is maintained 
In the far field. The coupling of a continuous vortex-sheet representa- 
tion and a concentrated vortex-line representation for solving the non- 
linear lifting-surface problem is called the "Nonlinear Hybrid-Vortex 
(NHV)* method. 

STEADY-FLOW PROBLEM-Equation (12) requires the flow to be tangent 
to w while equation (13) requires this tangential flow to be parallel 
to the vorticity direction. Therefore, if the flow direction is forced 
to be parallel to the vorticity direction on the surface w, the boundary 
conditions of equations (12) and (13) are automatically satisfied. 

Next, we outline the method of solution. 

Once the wing and its free-shear layers are represented by vortex 
sheets, equations (10) and (15) are automatically satisfied. The basic 
unknowns in the present problem are the vorticity distribution ij and 
the free-vortex sheet w. They are determined by satisfying the remaining 
boundary conditions, equations (H)-(14), through a finite-element type 
approach. 

In this approach, the bound-vortex sheet (representing the wing) is 
divided into quadrilateral bound-vortex panels while the free-vortex 
sheets (representing the free-shear layers) are divided into triangular 
free-vortex panels, see Figure 30. On each vortex panel (bound or free), 
a local vorticity distribution with undetermined coefficients is defined 
In a local -coordinate system (Figures 31, 32). The local distribution 
Is selected such that the solenoidal property of vorticity is satisfied. 
The continuity of vorticity (a compatibility condition) is enforced at 
certain nodes on the interelement boundaries of adjacent panels. 

The remaining boundary conditions, equations (11)-(14), are enforced 
at certain nodes of the vortex panels to obtain the undetermined coef- 
ficients of the local vorticity distribution and the shape of the 
free-vortex panels. Kutta condition, equation (14), is enforced at the 
nodes of the bound-vortex panels along the edges of separation. The no- 
penetration condition, equation (11), is enforced at the average points 
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of the bound-vortex panels. The no-penetration and no-pressure- jump 
conditions* equations (12) and (13), are simultaneously satisfied at 
the nodes of the free-vortex panels, Figure 33. 

To satisfy these conditions, an Iterative technique is followed 
which alternatively yields the local vorticity distribution on the 
bound-vortex panels and the shape of the free-vortex panels. Figure 34. 
During a typical iterative cycle, an overdetermined set of algebraic 
equations are solved for the undetermined coefficients of the local 
vorticity distributions. This is followed by adjusting the surface w 
such that u and (e«, - V(j») at the nodes of the free-vortex panels become 
parallel. The overdetermined set of equations consists of the equations 
obtained from the continuity of vorticity condition, Kutta condition, 
the no-pen'‘'^'* 3 t;.';'i condition on the wing, and a symmetry condition 
applied at ihe noiics along the line of symmetry. This set is solved by 
a least-square technioue. 

Once the iterative technique converges, the pressure distribution 
Is calculated and this is followed by calculating the total aerodynamic 
characten?:tics. Convergence of the technique is expedited by using an 
Initial guess for v; provided by NOV-method. More details of the method 
are given in reference 12. 

UNSTEADY-FLOW PROBLEM- The source of flow undsteadiness can be a 
time dependent angle of attack or a time-dependent freestream speed. 

In this time-domain approach the function of angle of attack is divided 
Into discrete changes in the angle of attack corresponding to discrete 
changes in time; i.e., at t*to, a*ao and at t=to + At, ao+4a etc. The 
problem is then solved at each time step where the solution of each step 
serves as the initial condition for the next discrete time step. The 
case of an impulsively started wing from rest is considered for the 
purpose of explaining the details of the technique. In this case, we 
set a*0 and replace e» by -e* in equations (16), (17) and (19)-(22). 

The following steps explain the procedure to enforce the boundary condi- 
tions of equations (19)-(22) to obtain oi(r,t) and w(r,t): 
a. The initial condition, at t>tQ, is considered to be a wing without 
a wake surface having a velocity of -e«. At this instant, we assume 
that not enough time has passed for the vorticity to be convected from 
the separation edges. The initial vorticity distribution on the wing 
io(r,t) is obtained from the least-square solution of an over-determined 
set of linear algebraic equations in the unknown coefficients describing 
the vorticity distribution. The set of equations consists of the no- 
penetration conditions, the continuity of vorticity conditions, the 
Kutta conditions, and the symmetry conditions. 

It should be noted that the initial vorticity distribution on the 
wing is such that the circulation around any cloied curve embracing a 
wing section or equivalently the outflow of vorticity from a surface 
enclosed by this curve is zero. Consequently, a starting vortex of 
opposite strength to that of the vorticity on the wing develops at the 
edges of separation. Next, the starting vortex is convected downstream 
with the local particle velocity, the flow existing at the preceding 
Instant is disturbed, and the vorticity distribution on the wing changes 
creating a new starting vortex to be shed downstream. 

This process continues and by the end of the first time-step, at 
t*tQ+At, a free-vortex strip, consisting of triangular free-vortex panels 
attached to the separation edges, is created. The free-vortex strip 
obeys the conditions given by equations (20) and (21) and its upstream 
edge satisfies Kutta condition at the separation edges, equation (22). 


19 


b. At t«to+At, one needs v(r,to+At) to determine the width of the vortex 
strip. Since this velocity is unknown so far, a good estimate is taken 
as the velocity at the preceding time step v(r,to), which is completely 
known from the solution corresponding to the initial condition. Within 
the time step At, the displacement Ar of any panel node is found from 

A? - r(t^ + At) - rCt^j - At • V(r,tjj) (32) 

Equation (22) is the same as equation (28) of the NDV-method. 

At t»tg+At, one also needs the vorticity of the triangular vortex 
panels forming the free-vortex strip, iI>yy(r,to+At) . For each triangular 
panel, we express its five unknown coefficients, describing its linear 
vorticity distribution, in terms of the five unknown coefficients of 
the adjacent bound-vortex panel, at the separation edges, at t»to 
(the coefficients of the bound panel are already known) and at t»to+At 
(the coefficients of the bound panel are still unknown). This is 
achieved by satisfying the following conditions. At t»to+At the vorti- 
city is continuous at the global nodes (on the separation edge) between 
the bound- and free-vortex panels. Also, at t*to+At, the fluid particles 
along the downstream edges of the triangular panels are the same parti- 
cles which existed at the edge of separation at t=to. According to 
Helmholtz theorem [Du/Dt » (u>*^)v], the vorticity of these particles 
changes as they are convected downstream. But according to Kelvin 
theorem (Dr/Dt » 0), the circulation around these particles remains 
constant and hence additional equations are written to satisfy Kelvin 
theorem between t«tg and t=to+At. 

Next, the vorticity distribution in the wing w(f,tQ+At) is obtained 
from the overdetermined set of algebraic equations. 

c. At t»to+2At, a new free-vortex strip is created along the separation 
edges, the first shed free-vortex strip is convected downstream under 
the condition Dr/Dt * 0, and_step (b) is repeated to find locations of 
the free-panels nodes and al(r,to+2At) . 

d. The steady state is reached once the change in the vorticity distri- 
bution is less than a prescribed error. 

To calculate the net pressure coefficient at any point on the wing 
surface r at any time t^ » to + tj^ • At, we apply the unsteady Bernoulli's 
equation 

- . , . . 3A«(r,tJ 

ACp(r,tk) » -2Vj(r,t^) • [7^(r,t^) + (e^ + otxr)] - 2— 

(33) 

where Vj(r,tk) is the jump in the tangential velocity across the wing, 
which is completely known from the vorticity vector at this location, 
u(r,tk). 

NUMERICAL EXAMPLE FOR A STEADY-FLOW- The developed computer program 
of Figure 34 is used to solve for the steady flow past a rectangular wing 
having side-edge separation. The wing is of aspect ratio of one at 9.7° 
angle of attack. The wing is divided into 6x6 quadrilateral bound- 
vortex panels, the trailing edge free-vortex sheet is divided into 6x6 
quadrilateral free-vortex panels, and the side edge free-vortex sheet is 
divided into a 6 vortex strips and each is divided into a different 
number of quadrilateral free-vortex panels such that the last panel in 
each strip occupy the same chordwise station as that of the last panel 
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of the trail ing-edge free-vortex panel. 

Figure 35 shows a typical converged solution of the spanwise and 
chordwise components of vorticity at the local nodes of the quadrila- 
teral vortex panels. At any node, the upper number is the spanwise 
component ur and the lower number is the chordwise component ur. It 
is seen that the continuity of vorticity condition, Kutta condition 
and the symnetry condition are satisfied at the common nodes, at the 
nodes of the trailing and side edges, and at the notes of the line of 
symmetry; respectively. 

Figure 36 shows the spanwise variation of the section normal- 
force coefficient at three iteration steps. The solution converges in 
the third itera^on step and is in good agreement with the experimental 
data of Scholtz^^. The figure also shows the solution of the same case 
oLi:a1ned by the NDV-method with 6x6 bound-vortex lines. One can con- 
clude that the NDV-method underestimates the normal-force coefficient 
near the wing tip. If the number of bound-vortex lines of the iJDV- 
method is increased to 9x7 (Figure 16), the solution agrees with that 
obtained by the NHV-method with 6x6 bound-vortex panels. This clearly 
shows that a less number of vortex panels gives the same accuracy as 
that obtained by a large number of bound-vortex lines. 

Figure 37 shows the converged solution of the spanwise variation 
of the section pitching-moment coefficient for the wing. The results 
of the NDV-method with 6x6 and 9x7 bound-vortex lines are also included 
in the figure. 

Figure 38 shows the spanwise variation of the net-pressure coeffi- 
cient at different chord stations with different number of panels for 
a rectangular wing at 15® angle of attack. The corresponding results 
of the NDV-method are also shown in the figure. 

The present computer program is computationally efficient when it 
is compared with other existing codes which use high-order doublet 
distribution. The CPU time on the CYBER 175 for this case is about 
200 seconds. 

NUMERICAL EXAMPLE FOR AN UNSTEADY FLOW-The same computer program 
is also used to solve for the unsteady flow past impulsively strated 
wings from rest without side edge separation. In the two examples, 
the wing is divided into 5x5 quadrilateral bound-vortex panels with a 
sine distribution in the chordwise direction and a cosine distribution 
in the spanwise direction. In the present cases, the dimensionless 
time step is equivalent to 0.48 wh le the root-chord length is 5 units. 

Figure 39 shows the distribution of the lift coefficient for a 
rectangular wing of aspect ratio of three at 5® angle of attack for t»2 
and t*22. The present results are compared with the steady numerical 
data of reference (83) where 196 panels of constant potential function 
are used. It is also compared with the experimental data of Reference 
(82). Although we used 25 panels in the present case, the results 
compares well with the given data. 

Figure 40 compares the growth of indicial lift for the same wing 
considered above with the numerical data of reference (84) where 100 
panels of constant doublets are used. 

Figure 41 shows the distribution of net pressure coefficient for 
a rectangular wing of aspect ratio of two at 20.5® angle of attack for 
t-2, 10, 21 and 22. 

On the CYBER 175 computer, the CPU time for each case with 5x5 


bound panels and 22 time steps Is about 10 minutes. 

Currently, work Is underway to Increase the number of panels, 
use the efficient far field calculations, and calculate cases with 
side- and leading-edge separations. 

CONCLUDING REMARKS 

Steady and unsteady. Incompressible free-wake analysis has been 
presented in detail. In the first part of the lecture, the motive behind 
the problem is established. Some of the existing experimental measure- 
ments and data are discussed and the approximations for the mathematical 
model are deduced. Using the model, the general formulation of the 
problem and special steady and unsteady cases are developed. The 
existing methods of solution are divided into four main groups; 
slender-body methods, nonlinear vortex lattice methods, doublet-panel 
methods, and nonlinear vortex-panel methods. 

Two specific methods (one from the second group; the NDV -method 
and the other from the fourth group; the NHV-method) are presented In 
detail with numerical examples showing their capabilities and success. 

I would like to emp‘'=»sl 2 e that the nonlinear discrete-vortex method should 
not be abandoned In of the success of the MNDV-method. It Is simple, 
accurate, and efficient. For more accuracy, all we need Is to replace 
the near-field calculations by continuous vortex panels Instead of the 
vortex lines. 

The solution of the free-wake problem is far from being complete. 

The problem of ccxnpresslbll ity effects in high subsonic and transonic 
flow needs considerable efforts in the future. In the transonic regime, 
the problems of shock capturing and shock-free wake interaction will 
require careful analysis and bigger computing machines. Currently, 
several researchers are approaching the present problem using the 
finite difference method. It remains to be seen how successful this 
approach is. 
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Figure 1. Examples of Vortex Flows 
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Figure 2. Types of separation near the leading edge 
of highly swept back wings. 
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Figure 3. Lift Coefficient VZ angle of attack for 
highly swept back wings. 
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Figure 4. Inviscid model of primary separation 
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Figure 13. Typical solution of the free-vortex sheets for a canard- 
wing configuration, steady flow, AR*1,66, a=15", 5x5 and 
8x8 lattices for the can«,rd and main wing; respectively, 
NOV -met hod. 



Figure 14. Typical solution of the leading-edge vortex sheet for a 
wing-body combination (no leeward separation), AR=2, 
f*9.86, 0=12.15®, 9x9 lattice for wing, 25x16 lattice for 
the body. 
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Typical solution of the unsteady wake for a delta wing, 
AR-1, o-lS^+A sin t, 6x6 lattice, NDV-method. 
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Figure 18. Upper and lower, spanwise pressure variation on a steadily 
rolling delta wing, AR-0.7, a=0°, n »-0.2, 8x8 lattice, 
x/c^»0.778. ^ 
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Figure 19. Lift and pitching-moment coefficients for a wing-body 
combination, AR=2, f=9.86, 9x9 lattice for wing, 

25x16 lattice for body. 
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Figure 20. Variation of the normal-force and pitching-moment 

coefficients with increasing and^decreasing angle of 
attack, rectangular wing, AR*1, a*l, 4x4 lattice, 

NOV -method. 
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Fi9ur6 21. Viriation of the pi tchlnQ-moment and normal-force 

coefficients for a delta wing in pitching oscillation, 
AR-1, a-15"+4 sin w t. 4x4 lattice, NDV-method. 
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Figurt S5 . Chordwlst Variation of the 
Lift Coefficient at Various 
Spanwisc Station; for a 
Rectangular t^ing, AR-3, a*5* 

SxS Bound Panels, zero thickness. 
No Tip Separation. 



Figured. GrtMxth of Indicia! Lift for a • 
Rectangular King, ARa3, a*5*. 
SxS Bound Panels, zero Thick- 
ness. No Tip Separation. 



pressure coefficient at 
different chordwisc station 
for a Rectangular Wing, 
AR*2, a"20.S*, 5x5 Bound 
Panels, zero-ThIckness, No 
Tip Separation. 




